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Abstract 

In this review we present the main features of the current status of 
neutrino physics. After a review of the theory of neutrino mixing and 
oscillations, we discuss the current status of solar and atmospheric neu- 
trino oscillation experiments. We show that the current data can be 
nicely accommodated in the framework of three-neutrino mixing. We 
discuss also the problem of the determination of the absolute neutrino 
mass scale through Tritium /3-decay experiments and astrophysical ob- 
servations, and the exploration of the Majorana nature of massive neu- 
trinos through neutrinoless double-/? decay experiments. Finally, future 
prospects are briefly discussed. 
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1 Introduction 



The last five years have seen enormous progress in our knowledge of neutrino physics. 
We have now strong experimental evidences of the existence of neutrino oscillations, 
predicted by Pontecorvo in the late 50's [UE], which occur if neutrinos are massive and 
mixed particles. 

In 1998 the Super-Kamiokande experiment [Sj provided a model-independent proof 
of the oscillations of atmospheric muon neutrinos, which were discovered in 1988 by the 
Kamiokande |4, and 1MB [3] experiments. The values of the neutrino mixing parameters 
that generate atmospheric neutrino oscillations have been confirmed at the end of 2002 
by the first results of the K2K long-baseline accelerator experiment [Hj, which observed 
a disappearance of muon neutrinos due to oscillations. 

In 2001 the combined results of the SNO [Zj and Super-Kamiokande jSj experiments 
gave a model-independent indication of the oscillations of solar electron neutrinos, which 
were discovered in the late 60's by the Homestake experiment j^j. In 2002 the SNO 
experiment jTH] measured the total flux of active neutrinos from the sun, providing a 
model-independent evidence of oscillations of electron neutrinos into other flavors, which 
has been confirmed with higher precision by recent data ^1] . The values of the neutrino 
mixing parameters indicated by solar neutrino data have been confirmed at the end of 
2002 by the KamLAND very- long-baseline reactor experiment [12], which have measured 
a disappearance of electron antineutrinos due to oscillations. 

In this paper we review the currently favored scenario of three-neutrino mixing, which 
is based on the above mentioned evidences of neutrino oscillations. In SectionEJwe review 
the theory of neutrino masses and mixing, showing that it is likely that massive neutrinos 
are Majorana particles. In Section El we review the theory of neutrino oscillations in 
vacuum and in matter. In Section H] we review the main results of neutrino oscillation 
experiments. In Section E] we discuss the main aspects of the phenomenology of three- 
neutrino mixing, including neutrino oscillations, experiments on the measurement of the 
absolute scale of neutrino masses and neutrinoless double-/? decay experiments searching 
for an evidence of the Majorana nature of massive neutrinos. In Section E] we discuss 
some future prospects and in Section [7| we draw our conclusions. 

For further information on the physics of massive neutrinos, see the books in Refs. |T3*1 

□HUE], the reviews in Refs. [IEI El UHl HHI EIB EH E21 ESI EH ESI EEl EH EEl EHl EHl IHEIl 

EH E21 EH1 EH E3 EI3 EH EE] and the references in Ref . |H1 . 

2 Neutrino masses and mixing 

The Standard Model was formulated in the 60's [1011111112] on the basis of the knowledge 
available at that time on the existing elementary particles and their properties. In partic- 
ular, neutrinos were though to be massless following the so-called two-component theory 
of Landau [13], Lee and Yang and Salam [IS], in which the massless neutrinos are 
described by left-handed Weyl spinors. This description has been reproduced in the Stan- 
dard Model of Glashow [10], Weinberg [H] and Salam |12] assuming the non existence 
of right-handed neutrino fields, which are necessary in order to generate Dirac neutrino 
masses with the same Higgs mechanism that generates the Dirac masses of quarks and 
charged leptons. However, as will be discussed in Section HJ in recent years neutrino 
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experiments have shown convincing evidences of the existence of neutrino oscillations, 
which is a consequence of neutrino masses and mixing. Hence, it is time to revise the 
Standard Model in order to take into account neutrino masses (notice that the Standard 
Model has already been revised in the early 70's with the inclusion first of the charmed 
quark and after of the third generation). 

2.1 Dirac mass term 

Considering for simplicity only one neutrino field u, the standard Higgs mechanism gen- 
erates the Dirac mass term 

£ D = -m D vv = -m D (V^u L + T>Zv R ) , (2.1) 

with iiid = yv/y/2 where y is a dimensionless Yukawa coupling coefficient and v/y2 is 
the Vacuum Expectation Value of the Higgs field, vl and vr are, respectively, the chiral 
left-handed and right-handed components of the neutrino field, obtained by acting on v 
with the corresponding projection operator: 

1 - 7 5 1 + 7 5 

V = V L + Vr , V L = P L V , Vr = PrV , P L = — - — , P R = — - — , (2.2) 

such that P L P R = P R P L = 0, P| = P L , P 2 R = P R , since (7 s ) 2 = 1. Therefore, we have 

P L u L = u L , P L v R = 0, Prv l = 0, P r vr = v r . (2.3) 

It can be shown that the chiral spinors and vr have only two independent components 
each, leading to the correct number of four for the independent components of the spinor 
v. 

Unfortunately, the generation of Dirac neutrino masses through the standard Higgs 
mechanism is not able to explain naturally why the neutrino are more than five order 
of magnitude lighter than the electron, which is the lightest of the other elementary 
particles (as discussed in Section HJ the neutrino masses are experimentally constrained 
below about 1-2 eV). In other words, there is no explanation of why the neutrino Yukawa 
coupling coefficients are more than five order of magnitude smaller than the Yukawa 
coupling coefficients of quarks and charged leptons. 

2.2 Majorana mass term 

In 1937 Majorana [IHj discovered that a massive neutral fermion as a neutrino can be 
described by a spinor ip with only two independent components imposing the so-called 
Majorana condition 

^ = r, (2.4) 

where ip c = Cip = C7 ip* is the operation of charge conjugation, with the charge- 
conjugation matrix C defined by the relations Cj^C^ 1 = — 7 M , & = C _1 , C T = —C. 
Since Crf T C~ l = 7 s and 7 5 7^ + 7^7 5 = 0, we have 

p L ^ = o, p L r R = ii>R, p R r L = r L , p R r R = o. (2.5) 
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In other words, ip c L is right-handed and ip R is left-handed. 

Decomposing the Majorana condition (|2.4jl into left-handed and right-handed com- 
ponents, ipL + ipjt = if>1 + if> R , and acting on both members of the equation with the 
right-handed projector operator P R , we obtain 

*pR = r L - (2.6) 

Thus, the right-handed component ip R of the Majorana neutrino field if) is not indepen- 
dent, but obtained from the left-handed component if>i through charge conjugation and 
the Majorana field can be written as 

iP = iP L + iPl. (2.7) 



This field depends only on the two independent components of ipL. Using the constraint 
()2.6|) in the mass term (j2.1|) . we obtain the Majorana mass term 

C M = -^m M (^lip L + ^hiPl) , (2-8) 

where we have inserted a factor 1/2 in order to avoid double counting in the Euler- 
Lagrange derivation of the equation for the Majorana neutrino field. 



2.3 Dirac-Majorana mass term 

In general, if both the chiral left-handed and right-handed fields exist and are indepen- 
dent, in addition to the Dirac mass term (j2.1|) also the Majorana mass terms for and 
u R are allowed: 

^L=~lrn L (piu L + VEul) , = ~ m R (W R v R + Vru° r ) . (2.9) 

The total Dirac -l-Majorana mass term 



can be written as 



^ D+M = ~H ^)( mL mD )( UL c )+ H.c. • (2.11) 
2 v ' \m D m R J \v R ) 

It is clear that the chiral fields vl and v R do not have a definite mass, since they are 
coupled by the Dirac mass term. In order to find the fields with definite masses it is 
necessary to diagonalize the mass matrix in Eq. (j2.11j) . For this task, it is convenient to 
write the Dirac+Majorana mass term in the matrix form 

C D+M = - W L MN L + H.c. , (2.12) 



with the matrices 

M=( , . 4 , 

m D m R I \u R 



M=\ mL mD , N L =l VL e ) . (2.13) 



5 



The column matrix N L is left-handed, because it contains left-handed fields. Let us write 
it as 

N L = Un L , with n L =(^j, (2.14) 

where U is the unitary mixing matrix (W = U~ v ) and til is the column matrix of the 
left-handed components of the massive neutrino fields. The Dirac+Majorana mass term 
is diagonalized requiring that 

U T MU=( n ! 1 °V (2.15) 

with rrik real and positive for k — 1, 2. 

Let us consider the simplest case of a real mass matrix M. Since the values of itll 
and m R can be chosen real and positive by an appropriate choice of phase of the chiral 
fields v L and u R , the mass matrix M is real if to d is real. In this case, the mixing matrix 
U can be written as 

U = Op, (2.16) 
where O is an orthogonal matrix and p is a diagonal matrix of phases: 

„ / cos-# siniA /ft 0\ , n , _, 

=(-B in tf COBtfj • "=(0 J ' ' 217 » 

with \pk\ 2 = 1- The orthogonal matrix O is chosen in order to have 

T MO=( m ' 1 °,V (2.18) 







leading to 



2m D , 1 



tan 2§ = , to: 



m R -m L 7 "~ 2 * 1 2 



m L + m R ± \J (m L - m R f + 4to^ 



(2.19) 



Having chosen to l and to r positive, to' 2 is always positive, but to' x is negative if to^ > 
mi,m R . Since 

^M„=/ 'MO,= (« , (2.20) 

it is clear that the role of the phases pk is to make the masses rrik positive, as masses 
must be. Hence, we have p\ = 1 always, and p\ = 1 if m[ > or p\ = — 1 if m! x < 0. 
An important fact to be noticed is that the diagonalized Dirac+Majorana mass term, 

£ d+m = 1 ^ mk lfav kL + H.c. , (2.21) 

fc=l,2 

is a sum of Majorana mass terms for the massive Majorana neutrino fields 

v k = v kL + if kL (A; = 1,2). (2.22) 
Therefore, the two massive neutrinos are Majorana particles. 
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2.4 The see-saw mechanism 



It is possible to show that the Dirac+Majorana mass term leads to maximal mixing 
(6 = 7r/4) if rriL = rriR, or to so-called pseudo-Dirac neutrinos if m L and tur are much 
smaller that |77&d| (see Ref. However, the most plausible and interesting case is 

the so-called see-saw mechanism |47| 11*51 """"""] . which is obtained considering mi = and 
|771d| "C J71r. In this case 

(iTLd) 2 IJId 2 

mi ~ <C m D , m 2 ~ , tanw ~ <1, p 1 — —l. (2.23) 

m fi m^ 

What is interesting in Eq. ()2.23|) is that mi is much smaller than mD, being suppressed 
by the small ratio m^jmR. Since m 2 is of order ttir, a very heavy z/ 2 corresponds to a 
very light i/i, as in a see-saw. Since mo is a Dirac mass, presumably generated with the 
standard Higgs mechanism, its value is expected to be of the same order as the mass of a 
quark or the charged fermion in the same generation of the neutrino we are considering. 
Hence, the see-saw explains naturally the suppression of mi with respect to m D , providing 
the most plausible explanation of the smallness of neutrino masses. 

The smallness of the mixing angle •& in Eq. 1)2.23)1 implies that u±l — —vl and z/ 2 l — v R . 
This means that the neutrino participating to weak interactions practically coincides 
with the light neutrino Ui, whereas the heavy neutrino z/ 2 is practically decoupled from 
interactions with matter. 

Besides the smallness of the light neutrino mass, another important consequence of the 
see-saw mechanism is that massive neutrinos are Majorana particles, as we have shown 
above in the general case of a Dirac+Majorana mass term. This is a very important 
indication that strongly encourages the search for the Majorana nature of neutrinos, 
mainly performed through the search of neutrinoless double-/? decay. 

An important assumption necessary for the see-saw mechanism is mi = 0. Such 
assumption may seem arbitrary at first sight, but in fact it is not. Its plausibility follows 
from the fact that Ul belongs to a weak isodoublet of the Standard Model: 

L L = (£) . (2.24) 

Since Ul has third component of the weak isospin I3 = 1/2, the combination v c l vl = 
—vEC'ul in the Majorana mass term in Eq. (|2.9J) has I3 = 1 and belongs to a triplet. Since 
in the Standard Model there is no Higgs triplet that could couple to v c l vl in order to form 
a Lagrangian term invariant under a SU(2)i transformation of the Standard Model gauge 
group, a Majorana mass term for ul is forbidden. In other words, the gauge symmetries 
of the Standard Model imply mi = 0, as needed for the see-saw mechanism. On the 
other hand, uid is allowed in the Standard Model, because it is generated through the 
standard Higgs mechanism, and tur is also allowed, because vr and v r vr are singlets of 
the Standard Model gauge symmetries. Hence, quite unexpectedly, we have an extended 
Standard Model with massive neutrinos that are Majorana particles and in which the 
smallness of neutrino masses can be naturally explained through the see-saw mechanism. 

The only assumption which remains unexplained in this scenario is the heaviness of 
rriR with respect to mo- This assumption cannot be motivated in the framework of the 
Standard Model, because wir is only a parameter which could have any value. However, 
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there are rather strong arguments that lead us to believe that the Standard Model is a 
theory that describes the world only at low energies. In this case it is natural to expect 
that the mass uir is generated at ultra-high energy by the symmetry breaking of the 
theory beyond the Standard Model. Hence, it is plausible that the value of is many 
orders of magnitude larger than the scale of the electroweak symmetry breaking and of 
mo, as required for the working of the see-saw mechanism. 

2.5 Effective dimension- five operator 

If we consider the possibility of a theory beyond the Standard Model, another question 
regarding the neutrino masses arises: is it possible that a Lagrangian term exists at the 
high energy of the theory beyond the Standard Model which generates at low energy an 
effective Majorana mass term for viP. The answer is yes [213131112!: the operator with 
lowest dimension invariant 1 under SU(2)x x U(l)y that can generate a Majorana mass 
term for v L after electroweak symmetry breaking is the dimension-five operator 2 

(Ll a 2 $) (T 1 ($ T a 2 L L ) + H.c. , (2.25) 

where g is a dimensionless coupling coefficient and A4 is the high-energy scale at which 
the new theory breaks down to the Standard Model. The dimension-five operator in 
Eq. ()2.25|) does not belong to the Standard Model because it is not renormalizable. It 
must be considered as an effective operator which is the low-energy manifestation of the 
renormalizable new theory beyond the Standard Model, in analogy with the old non- 
renormalizable Fermi theory of weak interactions, which is the low-energy manifestation 
of the Standard Model. 

At the electroweak symmetry breaking 

rh f ( f )+ \ Symmetry Breaking / \ , , 

$= Uv ' wa) ' ( } 

the operator in Eq. ()2.25j) generates the Majorana mass term for ul in Eq. (|2.9J) . with 



m L = 




(2.27) 



This relation is very important, because it shows that the Majorana mass mi is suppressed 
with respect to v by the small ratio v/Ai. In other words, since the Dirac mass term mo 
is equal to v/y/2 times a Yukawa coupling coefficient, the relation ()2.27|) has a see-saw 
form. Therefore, the effect of the dimension-five effective operator in Eq. ()2.25|) does 
not spoil the natural suppression of the light neutrino mass provided by the see-saw 

1 Since the high-energy theory reduces to the Standard Model at low energies, its gauge symmetries 
must include the gauge symmetries of the Standard Model. 

2 In units where h = c = 1 scalar fields have dimension of energy, fermion fields have dimen- 
sion of (energy) 3 / 2 and all Lagrangian terms have dimension (energy) 4 . The "dimension- five" char- 
acter of the operator in Eq. (j2.25(l refers to the power of energy of the dimension of the operator 
{h\ (72 3>) C _1 ($ T (72 Ll), which is divided by the mass M. in order to obtain a Lagrangian term with 
correct dimension. 
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mechanism. Indeed, considering m^ ~ m^/m^ and taking into account that mi = | | , 
from Eq. (J2.19)) we obtain 

(2.28) 

Equations 1)2.27)1 and (j2.28j) show that the see-saw mechanism is operating even if is 
not zero, but it is generated by the dimension- five operator in Eq. (12.25)1 . On the other 
hand, if the chiral right-handed neutrino field vr does not exist, the standard see-saw 
mechanism cannot be implemented, but a Majorana neutrino mass can be generated 
by the dimension-five operator in Eq. ()2.25j) . and Eq. (J2.27)) shows that the suppression 
of the light neutrino mass is natural and of see-saw type. 

2.6 Three-neutrino mixing 

So far we have considered for simplicity only one neutrino, but it is well known from 
a large variety of experimental data that there are three neutrinos that participate to 
weak interactions: z/ e , v T . These neutrinos are called "active flavor neutrinos". From 
the precise measurement of the invisible width of the Z-boson produced by the decays 
Z — > v a v a we also know that the number of active flavor neutrinos is exactly three 

OL 

(see Ref. |S31), excluding the possibility of existence of additional heavy active flavor 
neutrinos 3 . The active flavor neutrinos take part in the charged- current (CC) and neutral 
current (NC) weak interaction Lagrangians 

2 (2.29) 

^ Va~LlpVaL-, (2.30) 
a=e,fi,r 

where and are, respectively, the charged and neutral leptonic currents, $w is the 
weak mixing angle (sin 2, # w ~ 0.23) and g = e/ sin^w (e is the positron electric charge). 

Let us consider three left-handed chiral fields u eL , u^ L , v tL that describe the three 
active flavor neutrinos and three corresponding right-handed chiral fields v Si r, v S2 Ri u s 3 r 
that describe three sterile neutrinos 4 , which do not take part in weak interactions. The 
corresponding Dirac+ Majorana mass term is given by Eq. (j2.1()j) with 

r^-^i^Mj^ + H.c, (2.31) 

s,/3 

£l = -\ E^ M i m + H - c - > ( 2 - 32 ) 

a,f3 
s,s' 

3 More precisely, what is excluded is the existence of additional active flavor neutrinos with mass 
< mz/2 ~ 46 GeV [HI]. For a recent discussion of the possible existence of heavier active flavor neutrinos 
see Ref. [55] , 

4 Let us remark, however, that the number of sterile neutrinos is not constrained by experimental 
data, because they cannot be detected, and could well be different from three. 



mi 



m L - 



{m D y 



m R 



cc 



J, CC ^ p + H.c, with 



2V2 Jp 

NC 9 



2 cost? 



w 



• Jp NC ^, with Jp NC 
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where M D is a complex matrix and M L , M R are symmetric complex matrices. The 
Dirac+Majorana mass term can be written as the one in Eq. (|2.12jl with the column 
matrix of left-handed fields 



V R 




N L =\J C ) , with u L and u% = I v c S2R I , (2.34) 



and the 6x6 mass matrix 

' M L (M D ^ T 



s 2 J 

■ V KR, 



M = [ M ° M« J ' < 2 ' 35 > 

The mass matrix is diagonalized by a unitary transformation analogous to the one in 
Eq. flni: 

N L =Yn L , with n L = : , (2.36) 

\vblJ 

where V is the unitary 6x6 mixing matrix and n^L are the left-handed components of 
the massive neutrino fields. The mixing matrix V is determined by the diagonalization 
relation 

Y T MY = diag(mi,...,m 6 ) , (2.37) 

with rrik real and positive for k — 1, . . . , 6 (see Ref. J7j for a proof that it can be done). 
After diagonalization the Dirac+Majorana mass term is written as 

I 6 

/: D+M = --^m fc ^ i z/ fcL + H.c, (2.38) 

k=l 

which is a sum of Majorana mass terms for the massive Major ana neutrino fields 

Vk = VkL + v c kL (A; = 1, ... ,6) . (2.39) 

Hence, as we have already seen in Section 12.31 in the case of one neutrino generation, a 
Dirac+Majorana mass term implies that massive neutrinos are Majorana particles. The 
mixing relation ()2.36|) can be written as 

6 6 

v a L = ^2 Vafc VkL ( a = e ' ^' r ) ' u sR = Vsfc UkL ( s = Sl ' S2 ' s ^ ' ( 2A{) ) 

k=l k=l 

which shows that active and sterile neutrinos are linear combinations of the same massive 
neutrino fields. This means that in general active-sterile oscillations are possible (see 
Section EJ. 

The most interesting possibility offered by the Dirac+Majorana mass term is the 
implementation of the see-saw mechanism for the explanation of the smallness of the 
light neutrino masses, which is however considerably more complicated than in the case 
of one generation discussed in Section 12.41 Let us assume that M L = 0, in compliance 



10 



with the gauge symmetries of the Standard Model and the absence of a Higgs triplet 5 . 
Let us further assume that the eigenvalues of M R are much larger than those of M D , as 
expected if the Majorana mass term ()2.33j) for the sterile neutrinos is generated at a very 
high energy scale characteristic of the theory beyond the Standard Model. In this case, 
we can write the mixing matrix V as 

V = WU, (2.41) 

where both W and U are unitary matrices, and use W for an approximate block-diagonalization 
of the mass matrix M at leading order in the expansion in powers of (M^^M 13 : 



W T MW~ f M " sht ° ) . (2.42) 
U IVl] 



heavy 



The matrix W is given by 



(M v y(M R (M R )Y 1 M D 2(M D )t(M /? )t~ 1 

-2{M R y l M D {M R )~ l M D {M D y{M R y~ 



1 



(2.43) 



and is unitary up to corrections of order (M R ) 1 M D . The two 3x3 mass matrices Mjight 
and M heavy are given by 



'light 



-(M D ) T (M i? )" 1 M D , M heavy ~M R . (2.44) 



Therefore, the see-saw mechanism is implemented by the suppression of Mu g ht with re- 
spect to M D by the small ratio (M D ) T (M R ) _1 . The light and heavy mass sectors are 
practically decoupled because of the smallness of the off-diagonal 3x3 block elements in 
Eq. (fUlI). 

For the low-energy phenomenology it is sufficient to consider only the light 3x3 mass 
matrix Mii g ht which is diagonalized by the 3x3 upper-left submatrix of U that we call 
U, such that 

U T M light U = diag(mi, m 2 , m 3 ) , (2.45) 

where mi, 7712,7713 are the three light neutrino mass eigenvalues. Neglecting the small 
mixing with the heavy sector, the effective mixing of the active flavor neutrinos relevant 
for the low-energy phenomenology is given by 



3 



v<xL = U ak v kL (a = e, fi, t) , (2.46) 

k=l 

where un, v 2 Li v zl are the left-handed components of the three light massive Majorana 
neutrino fields. This scenario, called "three-neutrino mixing", can accommodate the 
experimental evidences of neutrino oscillations in solar and atmospheric neutrino exper- 
iments reviewed in Section 0] The phenomenology of three-neutrino mixing is discussed 
in Sectional 



5 For the sake of simplicity we do not consider here the possible existence of effective dimension- five 
operators of the type discussed in Section 1231 which in any case do not spoil the effectiveness see-saw 
mechanism. 
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The 3x3 unitary mixing matrix U can be parameterized in terms of 3 2 = 9 parameters 
which can be divided in 3 mixing angles and 6 phases. However, only 3 phases are physical. 
This can be seen by considering the charged- current Lagrangian ()2.29j) 6 . which can be 
written as 

3 

v a=e,fi,r k=l 

in terms of the light massive neutrino fields v\~ (k = 1,2,3). Three of the six phases in 
U can be eliminated by rephasing the charged lepton fields e, /i, r, whose phases are 
arbitrary because all other terms in the Lagrangian are invariant under such change of 
phases (see Refs. jSHEIlEB] and the appendices of Refs. [201 EHj)- On the other hand, the 
phases of the Majorana massive neutrino fields cannot be changed, because the Majorana 
mass term in Eq. ()2.38|) are not invariant 7 under rephasing of v^. Therefore, the number 
of physical phases in the mixing matrix U is three and it can be shown that two of these 
phases can be factorized in a diagonal matrix of phases on the right of U. These two 
phases are usually called "Majorana phases", because they appear only if the massive 
neutrinos are Majorana particles (if the massive neutrinos are Dirac particles these two 
phases can be eliminated by rephasing the massive neutrino fields, since a Dirac mass 
term is invariant under rephasing of the fields). The third phase is usually called "Dirac 
phase", because it is present also if the massive neutrinos are Dirac particles, being the 
analogous of the phase in the quark mixing matrix. These complex phases in the mixing 
matrix generate violations of the CP symmetry (see Refs. [T3*l IT4*1 ITT| I21j). 
The most common parameterization of the mixing matrix is 

U = R 23 Wl3 R\2 D(\ 2 i, A31 

'1 \ / r ,;; .s i; , ( 

C23 s 23 
v — s 23 C23 

C12C13 S12C13 
-S12C23 - c 12 s 23 s 13 e lVl3 c 12 c 23 - s l2 s 23 s 13 e lLpVi 
S12S23 - ci 2 c 23 s 13 e t{p13 -c 12 s 23 - s 12 c 23 s 13 e lipi3 



with Cij = cos^jj, Sij = sin^jj, where $i 2 , $23, $13 are the three mixing angles, (p i3 is the 
Dirac phase, A21 and A31 are the Majorana phases. In Eq. ()2.48j) Rij is a real rotation in 
the i-j plane, W13 is a complex rotation in the 1-3 plane and -D(A2i, A31) is the diagonal 
matrix with the Majorana phases. 

Let us finally remark that, although in the case of Majorana neutrinos there is no 
difference between neutrinos and antineutrinos and one should only distinguish between 

6 Unitary mixing has no effect on the neutral-current weak interaction Lagrangian, which is diagonal 

3 

in the massive neutrino fields, C^ c = —— — — — ^""^ vul 1 P ^kh Z p (G1M mechanism). 

^ COS l/\V 





fe=l 



7 In Field Theory, Noether's theorem establishes that invariance of the Lagrangian under a global 
change of phase of the fields corresponds to the conservation of a quantum number: lepton number L for 
leptons and baryon number B for quarks. The non-invariance of the Majorana mass term in Eq. (|2.38|) 
under rephasing of VkL implies the violation of lepton number conservation. Indeed, a Majorana mass 
term induces \AL\ = 2 processes as neutrinoless double-/3 decay (see Refs. [LHIlHIT7IOT1I^T] L 
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states with positive and negative helicity, it is a common convention to call neutrino 
a particles created together with a positive charged lepton and having almost exactly 
negative helicity, and antineutrino a particles created together with a negative charged 
lepton and having almost exactly positive helicity. This convention follows from the 
fact that Dirac neutrinos are created together with a positive charged lepton and almost 
exactly negative helicity, and Dirac antineutrinos are created together with a negative 
charged lepton and almost exactly positive helicity. 



3 Theory of neutrino oscillations 

In order to derive neutrino oscillations it is useful to realize from the beginning that 
detectable neutrinos, relevant in oscillation experiments, are always ultrarelativistic par- 
ticles. Indeed, as discussed in Section 01 the neutrino masses are experimentally con- 
strained below about 1-2 eV, whereas only neutrinos more energetic than about 200 keV 
can be detected in: 

1. Charged current weak processes which have an energy threshold larger than some 
fraction of MeV. For example 8 : 

• E th = 0.233 MeV for u e + 71 Ga -> 71 Ge + e" in the GALLEX [61J, SAGE |S2| 
and GNO [HS] solar neutrino experiments. 

• E t h = 0.81 MeV for u e + 37 C1 — > 37 At + e~ in the Homestake [§] solar neutrino 
experiment. 

• E t h = 1.8 MeV for v e +p — > n+e + in reactor neutrino experiments (for example 
Bugey |nH, CHOOZ jHS] and KamLAND [12]). 

2. The elastic scattering process v + e~ — > v + e~, whose cross section is proportional 
to the neutrino energy (cr(E) ~ <JoE/m e , with o"o ~ 10 _44 cm 2 ). Therefore, an 
energy threshold of some MeV's is needed in order to have a signal above the 
background. For example, E^ ~ 5 MeV in the Super-Kamiokande jtiti) EZ] solar 
neutrino experiment. 

The comparison of the experimental limit on neutrino masses with the energy threshold 
in the processes of neutrino detection implies that detectable neutrinos are extremely 
relativistic. 



3.1 Neutrino oscillations in vacuum 

Active neutrinos are created and detected with a definite flavor in weak charged- current 
interactions described by the Lagrangian 1)2.29)1 . The state that describes an active neu- 

8 In a scattering process v + A — > B + C the Lorentz-invariant Mandelstam variable s = (p„ + pa) 2 — 
('Pb + Pc) 2 calculated for the initial state in the laboratory frame in which the target particle A is at 
rest is s = 2Eitia + vn\- The value of s calculated for the final state in the center-of-mass frame is given 
by s = (Eb + Ec) 2 > (rriB +mc) 2 ■ Confronting the two expressions for s we obtain the neutrino energy 

(m B +m c ) 2 m A 



threshold in the laboratory frame Eth = 



2m A 2 
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trino with flavor a created together with a charged lepton a + in a decay process of type 9 

A^B + a + + u a (3.1) 

is given by 10 

3 

K) oc K> K« + bcclO> J^ B , (3.2) 

k=l 

where J^~* B is the current describing the A — > B transition. Neglecting the effect of neu- 
trino masses in the production process, which is negligible for ultrarelativistic neutrinos, 
from Eqs. flZZSD and (gHH) it follows that 

(u k ,a + \f cc \0)J^ B ^U: k . (3.3) 

Therefore, the normalized state describing a neutrino with flavor a is 

3 



fc=i 

This state describes the neutrino at the production point at the production time. The 
state describing the neutrino at detection, after a time T at a distance L of propaga- 
tion in vacuum, is obtained by acting on \u a ) with the space-time translation operator 11 

exp (—iET + iP • Lj, where E and P are the energy and momentum operators, respec- 
tively. The resulting state is 

ML,T)> = J2K k e- iEkT+ipkL k> , (3.5) 
fc=i 

where Ej. and p k are, respectively, the energy and momentum 12 of the massive neutrino 
Uf., which are determined by the process in which the neutrino has been produced. Using 
the expression of \i>k) in terms of the flavor neutrino states obtained inverting Eq. ()2.46j) . 

Wk) = ^2 u P h we ° btain 

k(l,t))= Yi (ib u ^ e ~ iEkT+ipkLu A\^i ( 3 - e ) 

(3=e,fi,T \fc=l / 

9 This is the most common neutrino creation process. Other processes can be treated with the same 
method, leading to the same result l|3.4(l for the state describing a ultrarelativistic flavor neutrino. 

10 The flavor neutrino fields v a are not quantizable because they do not have a definite mass and are 
coupled by the mass term. Therefore, the state \v a ) is not a quantum of the field v a . It is an appropriate 
superposition of the massive states \vk), quanta of the respective fields v kl which describes a neutrino 
created in the process (I3.1f) |55| . 

n We consider for simplicity only one space dimension along neutrino propagation. 

12 Since the energy and momentum of the massive neutrino z/j, satisfy the relativistic dispersion relation 
Ef. = pf. + m 2 ,, elementary dimensional considerations imply that at first order in the contribution of 

2 2 

the mass mk we have E k — E + £ — | and pk — E — (1 — £) — |, where E is the neutrino energy in the 

2E '2 /'. 

massless limit and ^ is a dimensionless quantity that depends on the neutrino production process. 
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which shows that at detection the state describes a superposition of different neutrino 
flavors. The coefficient of \vp) is the amplitude of v a — > vp transitions, whose probability 
is given by 



(L,T) = \(vp\v a (L,T))\ 1 < 



E 

k=l 



-iE k T+i Pk L jj^ 



(3.7) 



The transition probability ()3.7)1 depends on the space and time of neutrino propagation, 
but in real experiments the propagation time is not measured. Therefore it is necessary 
to connect the propagation time to the propagation distance, in order to obtain an ex- 
pression for the transition probability depending only on the known distance between 
neutrino source and detector. This is not a problem for ultrarelativistic neutrinos whose 
propagation time T is equal to the distance L up to negligible corrections depending on 
the ratio of the neutrino mass and energy 13 , leading to the approximation 



E k t - p k x ~ (E k - p k ) L 



El 



mi 



E k +Pk 



m k T 

E k + Pk IE 



(3. 



where E is the neutrino energy in the massless limit. This approximation for the phase of 
the neutrino oscillation amplitude is very important, because it shows that the phase of 
ultrarelativistic neutrinos depends only on the ratio m\L/E and not on the specific values 
of E k and p k , which in general depend on the specific characteristics of the production 
process. The resulting oscillation probability is, therefore, valid in general, regardless of 
the production process. 

With the approximation (|3.8|) . the transition probability in space can be written as 



Err* -imlL/2E TT 
U ak e * 1 Up k 

k 

Y,\U a k\ 2 \Up k \ 2 + 2ReY,Kk Upk U aj U* pj exp ( 



2E 



(3.9) 



where Am^ = m|— m|. Equation ()3.9|) shows that the constants of nature that determine 
neutrino oscillations are the elements of the mixing matrix and the differences of the 
squares of the neutrino masses. Different experiments are characterized by different 
neutrino energy E and different source-detector distance L. 
In Eq. (13. 9 j) we have separated the constant term 



p 

1 fa 



E 



V, 



ak 



\u, 



0k\ 



(3.10) 



from the oscillating term which is produced by the interference of the contributions of 
the different massive neutrinos. If the energy E or the distance L are not known with 
sufficient precision, the oscillating term is averaged out and only the constant flavor- 
changing probability (|3.10J) is measurable. 

13 A 

rigorous derivation of the neutrino transition probability in space that justifies the T — L approx- 
imation requires a wave packet description (see Refs. |69l 1271 170| and references therein). 
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In the simplest case of two-neutrino mixing 14 between u a , vp and V\, i>2, there is only 
one squared-mass difference Am = Am^ = m\ — m\ and the mixing matrix can be 
parameterized 15 in terms of one mixing angle 



— sin $ cos •& 

The resulting transition probability between different flavors can be written as 

P^(L) = sin 2 2tf sin 2 . (3.12) 

This expression is historically very important, because the data of neutrino oscillation 
experiments have been always analyzed as a first approximation in the two-neutrino 
mixing framework using Eq. ()3.12j) . The two-neutrino transition probability can also be 
written as 

J (Am 2 /eV 2 ) a/km)\ , x 

P» a ^(L) = sin 2 2tf sin 2 ( 1.27^ (e/G^Y J ' (3 " 13) 

where we have used typical units of short-baseline accelerator experiments (see below). 
The same numerical factor applies if L is expressed in meters and E in MeV, which are 
typical units of short-baseline reactor experiments. 

The transition probability in Eq. ()3.13|) is useful in order to understand the classifi- 
cation of different types of neutrino experiments. Since neutrinos interact very weakly 
with matter, the event rate in neutrino experiments is low and often at the limit of the 
background. Therefore, flavor transitions are observable only if the transition probability 
is not too low, which means that it is necessary that 

^> 0.1-1. (3,4) 

Using this inequality we classify neutrino oscillation experiments according to the ratio 
L/E which establishes the range of Am 2 to which an experiment is sensitive: 

Short-baseline (SBL) experiments. In these experiments L/E < leV~ 2 . Since the 
source-detector distance in these experiment is not too large, the event rate is 
relatively high and oscillations can be detected for Am 2 L/4E > 0.1, leading a 
sensitivity to Am 2 > 0.1 eV 2 . There are two types of SBL experiments: reactor v e 
disappearance experiments with L ~ 10 m, E ~ 1 MeV as, for example, Bugey 64J; 
accelerator experiments with L < 1 km, E > 1 GeV, as, for example, CDHS [2T| 
(vn — > fix), CCFR [12] {i>n -> ^, Vix — > fe and v e — > u T ), CHORUS [IS] (f M — > v T 
and v e z/ T ), NOMAD 1A\ (^ u T and -> i/ e ), LSND JZ3 (^ -»• v e and 



A' 



KARMEN [in] (^ 



Long-baseline (LBL) and atmospheric experiments. In these experiments L/E < 
10 4 eV~ 2 . Since the source-detector distance is large, these are low-statistics ex- 
periments in which flavor transitions can be detected if Am 2 L/AE > 1, giving a 

1 This is a limiting case of three-neutrino mixing obtained if two mixing angles are negligible. 
15 Here we neglect a possible Majorana phase, which does not have any effect on oscillations (see the 
end of Section . 
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sensitivity to Am 2 > 10 4 eV 2 . There are two types of LBL experiments analogous 
to the two types of SBL experiments: reactor v e disappearance experiments with 
L ~ 1km, E ~ lMeV (CHOOZ [Hj and Palo Verde j7B)): accelerator ^ experi- 
ments with L < 10 3 km, E > 1 GeV (K2K for ^ ^ and ^ -> z/ e , MINOS 
for z/ M — » z/^ and — > u e , CNGS jSH] for z/ M — > z/ r ). Atmospheric experiments 
(Kamiokande [EU, 1MB [12! , Super-Kamiokande 0, Soudan-2 jHSJ, MACRO jH3J) 
detect neutrinos which travel a distance from about 20 km (downward-going) to 
about 12780 km (upward-going) and cover a wide energy spectrum, from about 100 
MeV to about 100 GeV (see Section IP). 

Very long-baseline (VLBL) and solar experiments. The only existing VLBL is the 
reactor v e disappearance experiment KamLAND [12; with L ~ 180 km, E ~ 3 MeV, 
yielding L/E ~ 3 x 10 5 eV~ 2 . Since the statistics is very low, the KamLAND exper- 
iment is sensitive to Am 2 > 3 x 10~ 5 eV 2 . A sensitivity to such low values of Am 2 
is very important in order to have an overlap with the sensitivity of solar neutrino 
experiments which extends over the very wide range 10- 8 eV~ 2 < Am 2 < lO^eV" 2 
because of matter effects (discussed below). Solar neutrino experiments (Homes- 
take 0, Kamiokande [HS], GALLEX jEI], GNO SAGE |52|, Super-Kamiokande 
jHHl EZI, SNO [3 Uni HI]) can a l so measure vacuum oscillations over the sun- 
earth distance L ~ 1.5 x 10 8 km, with a neutrino energy E ~ 1 MeV, yielding 
L/E ~ 10 12 eV" 2 and a sensitivity to Am 2 > 10~ 12 eV 2 . 

3.2 Neutrino oscillations in matter 

So far we have considered only neutrino oscillations in vacuum. In 1978 Wolfenstein j86j 
realized that when neutrinos propagate in matter oscillations are modified by the coherent 
interactions with the medium which produce effective potentials that are different for 
different neutrino flavors. 

Let us consider for simplicity 16 a flavor neutrino state with definite momentum p, 

K(p)) = J2 U «kWk(p))- (3-15) 

k 

The massive neutrino states \vk{p)) with momentum p are eigenstates of the vacuum 
Hamiltonian 7i : 

Wo K(p)) = E k \vk{p)) , with E k = \Jp 2 + m 2 k . (3.16) 
The total Hamiltonian in matter is 

n = n + n I , with n I \u a (p))=v a \u a {p)) } (3.17) 

where V a is the effective potential felt by the active flavor neutrino v a (a = e, /i, r) 
because of coherent interactions with the medium due to forward elastic weak CC and 

16 A more complicated wave packet treatment is necessary for the derivation of neutrino oscillations in 
matter taking into account different energies and momenta of the different massive neutrino components 
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z 




Figure 1: Feynman diagrams of the coherent forward elastic scattering processes that 
generate the CC potential Vcc through W exchange and the NC potential Vnc through 
Z exchange. 

NC scattering whose Feynman diagrams are shown in Fig. ^ The CC and NC potential 
are 



V cc = V2G F N e , V NC = ~\G F N n , (3.18) 

where G F is the Fermi constant, and N e and N n are, respectively, the electron and neutron 
number densities. As shown in Fig. ^ the CC potential Vcc is felt only by the electron 
neutrino, whereas the NC potential is felt equally by the three active flavor neutrinos. 
Moreover, since the NC potential due to scattering on electrons and protons are equal and 
opposite, they cancel each other (the medium is assumed to be electrically neutral) and 
only the NC potential due to scattering on neutrons contributes to Vnc- Summarizing, 
we can write 

V a = Vcc $ae + Vnc • (3.19) 

For antineutrinos the signs of all potentials are reversed. 

In the Schrodinger picture the neutrino state with initial flavor a obeys the evolution 
equation 

i — \v a (p,t)) = H\v a (p,t)), with \u a (p, 0)) = \p a {p)) ■ (3.20) 

Let us consider the amplitudes of v a — > v@ flavor transitions 

ipapfat) = (vp{p)\v a {p,t)) , with ip a p{p, 0) = 8 aP . (3.21) 

In other words, ip a /3(p,t) is the probability amplitude that a neutrino born at t — with 
flavor a is found to have flavor (3 after the time t. 

From Eqs. ()3.1fi|) . ()3.17j) and (J3.2()|) . the time evolution equation of the flavor transition 
amplitudes is 



E (E U ^ Ek U U + S Pp ^ P ip, t) . (3.22) 



Considering ultrarelativistic neutrinos for which 



2 

E k ~E + -±, p = E, t = x, (3.23) 
Ah/ 
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we have the evolution equation in space 
. d , , s / m? „ , \ , „ , / TT Am? 



i 



' P \ k 



MX) = (p+T^ + Vkc) MX) + V V fyfc -— ^ ^ + he he V CC ^ap(x) 



^rapy^) - Y 2E J ^ ' 7 J \ 4^ p 2E pk 

(3.24) 

where we put in evidence the term (p + m\/2E + VncOVw^ 3 -) which generates a phase 
common to all flavors. This phase is irrelevant for the flavor transitions and can be 
eliminated by the phase shift 

^Mx) -> ^(x) e^HH 1 -^"^' 1 '^ , (3.25) 
which does not have any effect on the probability of u a z/g transitions, 

P^>) = hM*)| 2 - (3.26) 
Therefore, the relevant evolution equation for the flavor transition amplitudes is 

1 d^ = S (E ^ ^# ^ + ^ J iM*) , (3.27) 

which shows that neutrino oscillation in matter, as neutrino oscillation in vacuum, de- 
pends on the differences of the squared neutrino masses, not on the absolute value of 
neutrino masses. Equation ()3.27j) can be written in matrix form as 

i-^ *a = ^ (U AM 2 C/t + A) V Q , (3.28) 
with, in the case of three-neutrino mixing, 



'0 

= ( I > Am2 = I Am 2i I > A = | | , (3.29) 

hi j 




,0 Am 2 





Mcc 







)• H 












V o 








where 

A cc = 2EV CC = 2V2EG F N e . (3.30) 

Since the case of three neutrino mixing is too complicated for an introductory discus- 
sion, let us consider the simplest case of two neutrino mixing between u e , and v\, z/ 2 . 
Neglecting an irrelevant common phase, the evolution equation ()3.28|) can be written as 

■d/^eA 1 /-Am 2 cos2tf + 2A CC Am 2 sin 2d\ f tfj ee \ 



dx \^ e /v \ Am 2 sin 2$ Am 2 cos 2$ J \ip efJ ,; 

where Am 2 = m 2 , — m\ and d is the mixing angle, such that 

v e = cos $ v\ + sin v-i , f M = — sin i? + cos i9 z/ 2 • (3.32) 

If the initial neutrino is a u e , as in solar neutrino experiments, the initial condition for 
the evolution equation (|3.3H) is 



i>ee(0)\ = (I 

Mo) J V° 
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(3.33) 



and the probabilities of u e — > 



Ufj, transitions and u e survival are 



1-P„ 



(X 



(3.34) 



In practice the evolution equation of the flavor transition amplitudes can always be 
solved numerically with sufficient degree of precision given enough computational power. 
Let us discuss the analytical solution of Eq. ()3.31|) in the case of a matter density profile 
which is sufficiently smooth. This solution is useful in order to understand the qualitative 
aspects of the problem. 

The effective Hamiltonian matrix in Eq. (j3.31|) can be diagonalized by the orthogonal 
transformation 



with ^ 



?Pee 



M 



cos $m sin t?M 
- sin t? M cos $m 



fa 

02 



(3.35) 

where ipk can be thought of as the amplitude of the effective massive neutrino in 
matter (although such probability is not measurable, because only flavor neutrinos can 
be detected). The angle $m is the effective mixing angle in matter, given by 



tan I'd 



tan I'd 



M 



.4 



cc 



Am 2 cos 2$ 



The interesting new phenomenon, discovered by Mikheev and Smirnov in 1985 
beautifully explained by Bethe in 1986 [HO]) is that there is a resonance for 

A C c = Am 2 cos 2$ , 

which corresponds to the electron number density 

Am 2 cos 2d 



N. 



2 V / 2£G I 



(3.36) 

(and 
(3.37) 

(3.38) 



In the resonance the effective mixing angle is equal to 45°, i.e. the mixing is maximal, 
leading to the possibility of total transitions between the two flavors if the resonance 
region is wide enough. This mechanism is called "MSW effect" in honor of Mikheev, 
Smirnov and Wolfenstein. 

The effective squared-mass difference in matter is 



Am 2 



M 



(Am 2 cos 2$ - A cc f + (Am 2 sin 2tf ) 2 



(3.39) 



Neglecting an irrelevant common phase, the evolution equation for the amplitudes of 
the effective massive neutrinos in matter is 



d fi/j 



dx \1p2 



1 

4E 



-Am 2 ! 





Am^ 



+ 



— id^fvi/dx 
id$M/dx 



01 

V>2 



with the initial condition 

V>i(o) 
Mo) 



cos t?^ — sin 



sin 0gj 



COS 



0? 



M 



sin ^ 



(3.40) 



(3.41) 
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Figure 2: Effective mixing angle (left) and effective squared-mass difference Am^ 
(right) in matter as functions of the electron number density N e divided by the Avogadro 
number N A , for Am 2 = 7 x l(T 6 eV 2 , sin 2 2$ = 10~ 3 . iV e R = Am 2 cos 2$/2V2EG F is the 
electron number density at the resonance, where $m = 45°. 



where i?^ is the effective mixing angle in matter at the point of neutrino production. 

If the matter density is constant, d$M/dx = and the evolutions of the amplitudes of 
the effective massive neutrinos in matter are decoupled, leading to the transition proba- 
bility 

P^(x) = sin 2 2tf M sin 2 (^^j , (3-42) 

which has the same structure as the two-neutrino transition probability in vacuum (|3.12j) . 
with the mixing angle and the squared-mass difference replaced by their effective values 
in matter. 

If the matter density is not constant, it is necessary to take into account the effect of 
d$ M /dx, 

d$ M 1 Am 2 sin 2$ dA cc 



dx 2 (Am 2 cos 2$ - A cc f + (Am 2 sin 2i?) 2 dx 
which is maximum at the resonance, 



dtf 



M 



dx 



1 d In N e 



R 2 tan 2$ dx 



(3.43) 



(3.44) 



This is illustrated in the left panel of Fig. Elfor Am 2 = 7 x 1CT 6 eV 2 , sin 2 2$ = 1(T 3 . One 
can see that for N e <C the effective mixing angle is practically equal to the mixing 
angle in vacuum, $m ~ i?, for N e ~ Nf" the effective mixing angle varies very rapidly 
with the electron number density, passing through 45° at N e = and going rapidly to 
90° for N e > Nf. 

The right panel of Fig. 121 shows the corresponding behavior of the effective squared- 
mass difference Amy, which is useful in order to understand how the presence of a 
resonance can induce an almost complete v e — > conversion of solar neutrinos. If the 
mixing parameters are such that at the center of the sun N e 3> A^^, the effective mixing 
angle is practically 90° and electron neutrinos are produced as almost pure z/ 2 . As the 
neutrino propagates out of the sun, it crosses the resonance at N e = N^, where the 
energy gap between V\ and is minimum. If the resonance is crossed adiabatically, the 
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neutrino remains v 2 and exits the sun as v 2 = sin u e + cos •& z/ M , which is almost equal to 
Ufj, if the mixing angle is small, leading to almost complete u e — > z/ M conversion. This is 
the case in which the MSW effect is most effective and striking, since a large conversion 
is achieved in spite of a small mixing angle. 

If the resonance is not crossed adiabatically, z/ 2 — ► V\ transitions occur in an interval 
around the resonance and the neutrino emerges out of the sun as a mixture of v 2 and u\, 
leading to partial conversion of u e into v^. Quantitatively, we can write the amplitudes 
of V\ and v 2 at any point x after resonance crossing as 



cos < exp U j ^ } da;' J A\ + sin < exp l-i J ^ } d,M .4* 



* Am?, 



^(z) 



R 

22 



n0 /. f XR Am 2 (x') A R . n0 / . PAm^x') , A , 
cos < exp li J ^ ! dx'j Af 2 + sin < exp l-i J ^ J dxM .4 

xexpf-z Am^jxO A ^ 



where 4.j^- is the amplitude of — > z/j transitions in the resonance. 

Considering x as the detection point on the earth, practically in vacuum, the proba- 
bility of u e survival is given by 

P Ue -*v e (x) = \ip ee (x)\ 2 , with ipee(x) = cos ipi(x) + sm$ijj 2 (x) . (3.47) 

If Am 2 > ICr 10 eV 2 all the phases in Eqs. and (f3~47?jl are very large and rapidly 

oscillating as functions of the neutrino energy. In this case, the average of the transition 
probability over the energy resolution of the detector washes out all interference terms 
and only the averaged survival probability 

PZ u e = cos 2 cos 2 < K| 2 + cos 2 sin 2 < |4*| 2 

+ sin 2 ^ cos 2 <|^f 2 | 2 + sin 2 ^sin 2 ^|^ R 2 | 2 , (3.48) 

which is independent from the sun-earth distance, is measurable. Taking into account 
that conservation of probability implies that 

I aR I 2 - I A R I 2 - 1 _ P I A R I 2 - I A* I 2 - P 

l^lll — 1^221 ~ 1 r c j 1^121 ~~ l* Al 2ll ~~ r ci {0.<±i!) 

where P c is the v\ v 2 crossing probability at the resonance, we obtain the so-called 
Parke formula [HI] for the averaged u e survival probability: 



P^ Ve = - + I - - P c j cos 2< cos 2tf . (3.50) 

This formula has been widely used for the analysis of solar neutrino data. 

The main problem in the application of the Parke formula ()3.50j) is the calculation of 
the crossing probability. This probability must involve the energy gap Am 2 A /2E between 
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V\ and v-i and the off diagonal terms proportional to d$M/dx in Eq. (|3.40j) . which cause 
the z/i ^ v 2 transitions. Indeed, the crossing probability can be written as |H2l ESI EH EH] 



exp (-f 7 F) - exp (-f 



I - exp (-f 7-^) 



(3.51) 



where 7 is the adiabaticity parameter 

Am 2 M /2E 



7 



2\d-d M /M 



Am 2 sin 2 2$ 



R 



2£cos 2^|dlniV e /dx| 



R 



If 7 is large, the resonance is crossed adiabatically and P c <C 1, leading to 



P 



sun, adiabatic 



1 

2 



^ + ^ cos 2$^ cos 2-d . 



(3.52) 



(3.53) 



The parameter F in Eq. ()3.51|) depends on the electron density profile. The left panel in 
Figure El shows the Standard Solar Model (SSM) electron density profile in the sun 
which is well approximated by the exponential 



N e (R) = N e (0)exp(-^- 



with JV e (0) = 245iV A /cm 3 , and R 







10.54 ' 
(3.54) 

where R is the distance from the center of the sun and R Q is the solar radius. For an 
exponential electron density profile the parameter F is given by [521 EH EH EH EH EHj 



1 - tan 2 . 



(3.55) 



For |d In iV e /dx| R the authors of Ref. |10()j suggested the practical prescription, verified 
with numerical solutions of the differential evolution equation, to calculate it numerically 
from the SSM electron density profile for R < O.9O4i? and take the constant value 
18.9/-Rq for R > O.9O4i? , where the exponential approximation ()3.54|) breaks down. 

For the analysis of solar neutrino data it is also necessary to take into account the 
matter effect along the propagation of neutrinos in the earth during the night (the so- 
called u v e regeneration in the earth"), which can generate a day-night asymmetry of the 
rates. The probability of solar v e survival after crossing the earth is given by [T%1 HUlj 



> sun+carth 



p: 



(l-SiCU) tet-sin 2 ^) 



cos2-# 



(3.56) 



Since the earth density profile is not a smooth function, the probability P^^ e must be 
calculated numerically. A good approximation is obtained by approximating the earth 
density profile with a step function (see Refs. |1U2| 11031 11U41 11051 Q06 ). According to 
Eq. ()3.40|) . the effective massive neutrinos propagate as plane waves in regions of constant 
density, with a phase exp (±zAm^Ax/4£?), where Ax is the width of the step. At the 
boundaries of steps the wave functions of flavor neutrinos are joined, according to the 
scheme 



(n) 



U M <$>( 



(n-1) 



23 




Figure 3: Left: Standard Solar Model electron density profile in the sun as a function of 
the ratio R/R Q 96J. The straight line represents the approximation in Eq. ()3.54jl . Right: 
The conventional names for regions in the tan 2 i5-Am 2 plane obtained from the analysis 
of solar neutrino data. The vertical dotted line correspond to maximal mixing. 



U M ^{x 2 -x l )Ul\ \UuHx 1 -x )Ul L ] U^(x ). (3.57) 



(2) 



(1) 



where xq is the coordinate of the point in which the neutrino enters the earth, x\, x 2 , ■ ■ ■ , 
x n are the boundaries of n steps with which the earth density profile is approximated, 
$(Ax) = diag(exp (— iAm^Ax /AE) , exp (iAm^Ax/AE)), and the notation [. . .L^ indi- 
cates that all the matter-dependent quantities in the square brackets must be evaluated 
with the matter density in the i th step, that extends from Xi-i to 2j. 

The right panel in Fig. El shows the conventional names for regions in the tan 2 d-Am 2 
plane obtained from the analysis of solar neutrino data. The Small Mixing Angle (SMA) 
region is the one where the mixing angle is very small and the resonant enhancement 
of flavor transitions due to the MSW effect is more efficient. However, as explained in 
Section l4~Tl there is currently a very strong evidence in favor of the Large Mixing Angle 
(LMA) region, in which both the mixing angle and Am 2 are large. Other regions with 
large mixing are: the low Am 2 (LOW) region, the Quasi- Vacuum-Oscillations (QVO) 
region, and the VACuum Oscillations region (VAC). In the SMA, LMA and LOW regions 
vacuum oscillations due to the sun-earth distance are not observable because the Am 2 is 
too high and interference effects are washed out by the average over the energy resolution 
of the detector (in these cases the Parke formula ([3. 50 J) applies). In the QVO region both 
matter effects and vacuum oscillations are important [107) 1108) 1109) 1100) . In the VAC 
region matter effects are negligible and vacuum oscillations are dominant. 

Concluding this Section on the theory of neutrino oscillations, let us mention that the 
evolution equation (J3.28)) allows to prove easily that the Majorana phases in the mixing 
matrix do not have any effect on neutrino oscillations in vacuum [110) llllj as well as in 
matter [112j . because the diagonal matrix of Majorana phases D(\2i, A31) on the right of 
the mixing matrix in Eq. ()2.48)) cancels in the product UAW?U^. Therefore, the Dirac 
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Table 1: Short-baseline experiments (SBL) whose data give the most stringent constraints 
on different oscillation channels. 

or Majorana nature of neutrinos cannot be distinguished in neutrino oscillations. 

4 Neutrino oscillation experiments 

In this Section we review the main results of the oscillation experiments which are con- 
nected with the existing model-independent evidences in favor of oscillations of solar and 
atmospheric neutrinos and the interpretation of the experimental data in the framework 
of three neutrino mixing, discussed in Sectional We do not discuss the results of several 
short-baseline neutrino (SBL) oscillation experiments, which have probed scales of Am 2 
bigger than about O.leV 2 , that are larger than the scales of Am 2 indicated by solar and 
atmospheric neutrino data. The SBL experiments whose data give the most stringent 
constraints on the different oscillation channels are listed in Table E 

All the SBL experiments in Table E did not observe any indication of neutrino oscil- 
lations, except the LSND experiment |114[ I75j. A large part of the region in the sin 2 2$- 
Am 2 plane allowed by LSND has been excluded by the results of other experiments 
which are sensitive to similar values of the neutrino oscillation parameters (KARMEN 
[TJj, CCFR fTT7j. NOMAD [JJ; see Ref. [HB! for an accurate combined analysis of LSND 
and KARMEN data). The MiniBooNE experiment [119] running at Fermilab will tell us 
the validity of the LSND indication in the near future. 

Some years ago the oscillations indicated by the LSND experiment could be accom- 
modated together with solar and atmospheric neutrino oscillations in the framework of 
four-neutrino mixing, in which there are three light active neutrinos and one light sterile 
neutrino (see Refs. [2~T ] I12()| I121| 1122] and references in Ref. [HE])- However, the global 
fit of recent data in terms of four- neutrino mixing is not good |123j . disfavoring such 
possibility. Therefore, in this review we discuss only three-neutrino mixing, which cannot 
explain the LSND indication, awaiting the response of MiniBooNE before engaging in 
wild speculations (see Refs. H2U EES EEE1 QUI EES EES CHOI ED]). 
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4.1 Solar neutrino experiments and KamLAND 

At the end of the 60's the radiochemical Homestake experiment jH] began the observation 
of solar neutrinos through the charged-current reaction |132t I133| 

v e + 37 C1 -> 37 Ar + e~ , (4.1) 

with a threshold = 0.814 MeV which allows to observe mainly 7 Be and 8 B neutrinos 
produced, respectively, in the reactions e~ + 7 Be — > 7 Li + v e (E = 0.8631 MeV) and 
8 B — > 8 Be* + e + + v e (E < 15 MeV) of the thermonuclear pp cycle that produces energy 
in the core of the sun (see Refs. |134t Ho]). 

The Homestake experiment is called "radiochemical" because the 37 Ar atoms were 
extracted every ~35 days from the detector tank containing 615 tons of tetrachloroethy- 
lene (C2CI4) through chemical methods and counted in small proportional counters which 
detect the Auger electron produced in the electron-capture of 37 Ar. As all solar neutrino 
detectors, the Homestake tank was located deep underground (1478 m) in order to have 
a good shielding from cosmic ray muons. The Homestake experiment detected solar elec- 
tron neutrinos for about 30 years measuring a flux which is about one third of the 
one predicted Standard Solar Model (SSM) 

(jjHom 

-1m =0.34 ±0.03. (4.2) 

This deficit was called "the solar neutrino problem". 

The solar neutrino problem was confirmed in the late 80 's by the real-time water 
Cherenkov Kamiokande experiment jHHj (3000 tons of water, 1000 m underground) which 
observed solar neutrinos through the elastic scattering (ES) reaction 

v + e~ -»• v + e" , (4.3) 

which is mainly sensitive to electron neutrinos, whose cross section is about six time larger 
than the cross section of muon and tau neutrinos. The experiment is called "real-time" 
because the Cherenkov light produced in water by the recoil electron in the reaction 
(14.3)1 is observed in real time. The solar neutrino signal is separated statistically from 
the background using the fact that the recoil electron preserves the directionality of the 
incoming neutrino. The energy threshold of the Kamiokande experiment was 6.75 MeV, 
allowing only the detection of 8 B neutrinos. After 1995 the Kamiokande experiment 
has been replaced by the bigger Super-Kamiokande experiment [HI H23 E2] (50 ktons of 
water, 1000 m underground) which has measured with high accuracy the flux of solar 8 B 
neutrinos with an energy threshold of 4.75 MeV, obtaining jHE] 

111 = 0.465 ±0.015. (4.4) 

^ES 

In the early 90's the GALLEX 61J (30.3 tons of 71 Ga, 1400 m underground) and 
SAGE [62] (50 tons of 71 Ga, 2000 m underground) radiochemical experiments started the 
observation of solar electron neutrinos through the charged-current reaction |135j 



. e ± 71 Ga^ 71 Ge±e-, (4.5) 
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which has the low energy threshold of 0.233 MeV, that allows the detection of the so-called 
pp neutrinos produced in the main reaction p + p — > d + e + + v e [E < 0.42 MeV) of the 
pp cycle, besides the 7 Be, 8 B and other neutrinos. After 1997 the GALLEX experiment 
has been upgraded, changing its name to GNO jHSl- The combined results of the three 
Gallium experiments confirm the solar neutrino problem: 

<1>G;1 0.56 ±0.03. (4.6) 



^Ga 



Although it was difficult to doubt of the Standard Solar Model, which was well tested 
by helioseismological measurements (see Ref. |136j ). and it was difficult to explain the 
different suppression of solar z/ e 's observed in different experiments with astrophysical 
mechanisms, a definitive model-independent proof that the solar neutrino problem is due 
to neutrino physics was lacking until the real-time heavy-water Cherenkov detector SNO 
[3nniE] (1 kton of D 2 0, 2073 m underground) observed solar 8 B neutrinos through the 
charged- current (CC) reaction 

u e + d-^p + p + e~, (4.7) 

with E^cc = 8.2 MeV and the neutral-current (NC) reaction 

u + d^p + n + u, (4.8) 

with £SNO c = 2.2 MeV, besides the ES reaction (@~3J) with £fNO = 7.0 MeV. The obser- 
vation of solar neutrinos through the CC and NC reactions has provided the breakthrough 
for the definitive solution of the solar neutrino problem in favor of new neutrino physics. 
The charged-current reaction is very important because it allows to measure with high 
statistics the energy spectrum of solar u e 's. The neutral current reaction is extremely 
important for the measurement of the total flux of active v e , and u T , which interact 
with the same cross section. 

In June 2001 the combination of the first SNO CC data [7] and the high-precision 
Super-Kamiokande ES data [8 allowed to extract a model-independent indication of the 
oscillations of solar electron neutrinos into active v^s and/or z/ T 's [7] (see also Refs. |137t 
138J). In April 2002 the observation of solar neutrinos through the NC and CC reactions 
allowed the SNO experiment ^0] to solve definitively the long-standing solar neutrino 
problem in favor of the existence of v e — > z/^, v T transitions. In this first phase [HI], called 
"D 2 phase", the neutron produced in the neutral-current reaction ()4.8j) was detected 
by observing the photon produced in the reaction 

n + d -> 3 H + 7 (£ 7 = 6.25 MeV) . (4.9) 

In September 2003 the SNO collaboration released the data obtained in the second phase 
[TT] . called "salt phase", in which 2 tons of salt has been added to the heavy water in 
the SNO detector, allowing the detection of the neutron produced in the neutral-current 
reaction (|4.8j) by observing the photons produced in the reaction 

n + 35 C1 -> 36 C1 + several 7 's (£* ot = 8.6 MeV) . (4.10) 

The better signature given by several photons and the higher cross-section of reaction 
(|4.10|) with respect to reaction f)4.9|) have allowed the SNO collaboration to measure with 



27 




10' 



10 



1 

tan 2 6 



10 



10 



" I I Rate excluded 
HMM1 Rate+Shape allowed 
M LMA 

Palo Verde excluded 
I Chooz excluded 



_l_ 



_l_ 



0.2 0.4 , 0.6 



0.8 



sin 



'20 



Figure 4: Left: Allowed regions of neutrino oscillation parameters obtained from the 
global analysis of solar neutrino data ^JJ. The best-fit point is marked by a star. Right: 
KamLAND excluded regions of neutrino oscillation parameters for the rate analysis and 
allowed regions for the combined rate and energy spectrum analysis at 95% CL [T2*] . 
At the top are the 95% CL. excluded region from CHOOZ |139j and Palo Verde [75] 
experiments, respectively. The dark area is the 95% CL. LMA allowed region obtained 
in Ref. 140J. The thick dot indicates the best fit of KamLAND data. 



good precision the total flux of active neutrinos coming from 8 B decay in the core of the 
sun [TT] : 

$sno = 5.21 ± 0.47 x 10 6 cm- 2 s _1 , (4.11) 
which is in good agreement with the value predicted by the Standard Solar Model (SSM) 

$SSM = g Q5 +1.01 x 1Q 6 cm -2 g -l (4 12) 

On the other hand, the flux of electron neutrinos coming from 8 B decay measured through 
the CC reaction K7\i is only [TT] 

$|NO = 1.59+0.10 x 1Q 6 cm -2 s -l (413) 

The fact that the ratio [TT] 

$sno 

-f^ = 0.306 ± 0.035 (4.14) 

differs from unity by about 19 standard deviations is a very convincing proof that solar 
electron neutrinos have transformed into muon and/or tau neutrinos on their way to the 
earth. 

The result of the global analysis of all solar neutrino data in terms of the sim- 
plest hypothesis of two-neutrino oscillations favors the so-called Large Mixing Angle 
(LMA) region with effective two-neutrino mixing parameters AmluN ~ 7 x 10 -5 eV 2 
and tan 2 $sun ~ 0.4, as shown in the left panel in Fig. HI taken from Ref. [TT] . 

A spectacular proof of the correctness of the LMA region has been obtained at the 
end of 2002 in the KamLAND long-baseline u e disappearance experiment i n which 
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the suppression 

7y KamLAND 

^°k"nd = 0-611 ±0-094- (4.15) 

expected 

of the v e flux produced by nuclear reactors at an average distance of about 180 km was 
observed. The right panel in Fig. 0] shows the regions of oscillation parameters allowed 
by KamLAND, compared with the allowed LMA region obtained in Ref. j!40| in 2002 
after the release of the data of the first D2O phase of the SNO experiment ^0]- From 
the right panel in Fig. |U one can see that the LMA region and the KamLAND allowed 
regions overlap in two subregions at Am§ UN ~ 7 x 10~ 5 eV 2 and Am^ UN ~ 1.5 x 10~ 4 eV 2 . 
Therefore, the combined fit of 2002 solar neutrino data and KamLAND data yielded two 
allowed LMA subregions. The 2003 SNO salt phase data lead to a restriction of the LMA 
region allowed by solar neutrino data, which favors the lower LMA subregion, as shown in 
the left panel in Fig. which depicts the most updated allowed region of the two-neutrino 
oscillation parameters obtained from the global analysis of solar and KamLAND neutrino 
data. The effective two-neutrino mixing parameters are constrained at 99.73% C.L. (3a) 
in the ranges |141j 

5.4 x 10" 5 eV 2 < Am| UN < 9.4 x 10~ 5 eV 2 , (4.16) 
0.30 < tan 2 $ S un < 0.64 , (4.17) 

with best-fit values 



Am 2 £ f N = 6.9 x 10~ 5 eV 2 , tan 2 tf^ N = 0.43 . (4.18) 

Maximal mixing is excluded at a confidence level equivalent to 5.4<r [TTj . 

Transitions of solar i/ e 's into sterile states are disfavored by the data. The right 
panel in Fig. El shows the allowed regions in the fB,totai^sin 2 r\ plane obtained in Ref. |142j 
before the release of the SNO salt data, where fB.totai = ^sb/^s™ is the ratio of the 
8 B solar neutrino flux and its value predicted by the Standard Solar Model (SSM) 
The parameter sin 2 ?? quantifies the fraction of solar i/ e 's that transform into sterile u s : 
v e — > cost] v a + sin 7] u s , where v a are active neutrinos. From the right panel of Fig. El it is 
clear that there is a correlation between fB.totai and sin 2 rj, which is due to the constraint 
on the total flux of 8 B active neutrinos given by the SNO neutral-current measurement: 
disappearance into sterile states is possible only if the 8 B solar neutrino flux is larger 
than the SSM prediction. The allowed ranges for $s B and sin 2 77 are [[142j 

$8 B = 1.00 ±0.06 $f| M , sin 2 77 < 0.52. (4.19) 

The allowed interval for <3>s B shows a remarkable agreement of the data with the SSM, 
independently from possible v e — > v s transitions. The recent SNO salt data do not allow 
to improve significantly the bound on sin 2 77 |143j . 

In the future it is expected that the KamLAND experiment will allow to reach a 
relatively high accuracy in the determination of Amg UN |144j . whereas new low-energy 
solar neutrino experiments or a new dedicated reactor neutrino experiment are needed 
in order to improve significantly our knowledge of the solar effective mixing angle $sun 
[Tl5llTl6irr47j . 
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Figure 5: Left: Allowed regions obtained from the global analysis of solar and KamLAND 
data P]. Right: Allowed 90%, 95%, 99%, 99.73% CL. regions obtained from the global 
analysis of solar and KamLAND data |142j . The best-fit points are marked by stars. 



4.2 Atmospheric neutrino experiments and K2K 

Atmospheric neutrinos are produced by cosmic rays (mainly protons) which interact with 
the atmosphere producing pions, which decay into muon and neutrinos, 



7T 



7T 



(4.20) 



At low energy the muons decay before hitting the ground into electrons and neutrinos, 



/i 1 



+ + V, 



ft 



+ v e + v u 



(4.21) 



Hence, the predicted ratio of z/ M + and v e + v e is about 2 at neutrino energy E < lGeV. 
At higher energies the ratio increases, but it can be calculated with reasonable accuracy 
(about 5%). On the other hand, the calculation of the absolute value of the atmospheric 
neutrino flux suffers from a large uncertainty (20% or 30%) due to the uncertainty of 
the absolute value of the cosmic ray flux and the uncertainties of the cross sections of 
cosmic ray interactions with the nuclei in the atmosphere (see Ref. |2Ej)- Therefore, 
the traditional way that has been followed for testing the atmospheric neutrino flux 
calculation is to measure the ratio of ratios 



[Nfa + u,)/N(u e + z/ e )] data 

Tat/ , — \ / A T / i — M 5 



R [N(u, + v,)/N{v e + u e )] theo 



(4.22) 



where the subscripts "data" and "theo" indicate, respectively, the measured and calcu- 
lated ratio. If nothing happens to neutrinos on their way to the detector the ratio of 
ratios should be equal to one. 

Atmospheric neutrinos are observed through high-energy charged-current interactions 
in which the flavor, direction and energy of the neutrino are strongly correlated with the 
measured flavor, direction and energy of the produced charged lepton. 

In 1988 the Kamiokande [3] and 1MB experiments measured a ratio of ratios 
significantly lower than one. The current values of R measured in the Super-Kamiokande 
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Figure 6: Left: Up-down Super-Kamiokande asymmetry as a function of momentum 
for e-like and /i-like events generated, respectively, by atmospheric u e , v e and 
[149 . The division of /i-like into fully contained (FC) and partially contained (PC) is 
explained, for example, in Ref . j2S] • The hatched region shows the theoretical expectation 
without neutrino oscillations. The dashed line for //-like events represents the fit of the 
data in the case of two-generation — ► v T oscillations with Am 2 = 3.5 x 10~ 3 eV 2 and 
sin 2 I'd = 1.0. Right: Allowed region contours for — > u T oscillations obtained by the 
Super-Kamiokande experiment |150j . 

experiment are [148 

R S ~ K (E < 1.33 GeV) = 0.638 ± 0.053 , (4.23) 
R S ~ K (E > 1.33 GeV) = 0.675 ± 0.087 . (4.24) 

The boundary of 1.33 GeV has been chosen by the Super-Kamiokande Collaboration for 
historical reasons connected with proton decay search. 

Also the Soudan-2 experiment |BB] observed a ratio of ratios significantly lower than 
one, 

i? Soudan - 2 = 0.69 ±0.12, (4.25) 

and the MACRO experiment [HI] measured a disappearance of upward-going muons. 

Although the values ()4.23j) . ()4.24|) and ()4.25j) of the ratio of ratios suggest an evidence 
of an atmospheric neutrino anomaly probably due to neutrino oscillations, they are not 
completely model- independent. 

The breakthrough in atmospheric neutrino research occurred in 1998, when the Super- 
Kamiokande Collaboration jH] discovered the up-down asymmetry of high-energy events 
generated by atmospheric u^s, providing a model independent proof of atmospheric 
disappearance. Indeed, on the basis of simple geometrical arguments the fluxes of upward- 
going and downward-going high-energy events generated by atmospheric z/^'s should be 
equal if nothing happens to neutrinos on their way from the production in the atmosphere 
to the detector (see Ref. |152j ). The last published value of the measured up-down 
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Figure 7: Left: 90% C.L. allowed region contours for — > v T oscillations obtained 
by the Super-Kamiokande, MACRO and Soudan-2 experiments - Right: Allowed 
region contours for disappearance obtained in the K2K experiment confronted with 
the allowed regions for — > v T oscillations obtained in the Super-Kamiokande experiment 

USB- 



asymmetry is [T 

= -0.31 ± 0.04 , (4.26) 

N% + Njfi™ J 

showing a 7a evidence of disappearance of atmospheric high-energy upward-going muon 
neutrinos. These neutrinos travel a distance from about 2650 to about 12780 km (0.2 < 
cos 8 < 1, where 9 is the nadir angle of the neutrino trajectory), whereas the downward- 
going neutrinos travel a distance from about 20 to about 100 km (—1 < cos# < —0.2). 
Therefore, the simplest explanation of the atmospheric neutrino data is neutrino oscil- 
lations. The left panel in Fig. El shows the Super-Kamiokande up-down asymmetry as 
a function of momentum for e-like and //-like events generated, respectively, by atmo- 
spheric u e , v e and v^. One can see that there is a clear deficit of high-energy upward 
going muon neutrinos with respect to downward-going ones, which do not have time to 
oscillate. On the other hand, there is no up-down asymmetry at low energies because 
also most of the downward-going muon neutrinos have time to oscillate and a possible 
asymmetry is washed out by a poor correlation between the directions of the incoming 
neutrino and the observed charged lepton (the average angle between the two directions 
is 55° at p = 400 MeV and 20° at 1.5 GeV 0). 

At the end of 2002 the long-baseline K2K experiment 6J confirmed the neutrino oscil- 
lation interpretation of the atmospheric neutrino anomaly observing the disappearance of 
accelerator z/^'s with average energy energy E ~ 1.3 GeV traveling 250 km from KEK to 
the Super-Kamiokande detector (only 56 of the 80.ll 5 ; 4 expected events were observed). 

The Super-Kamiokande atmospheric neutrino data and the data of the K2K exper- 
iment are well fitted by z/„ — * v T transitions with the effective two-neutrino mixing pa- 
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rameters constrained in the ranges |153j 

1.4 x l(T 3 eV 2 < Ami TM < 5.1 x 10~ 3 eV 2 , (4.27) 
sin 2 2$atm > 0.86 , (4.28) 

at 99.73% C.L. (3a), with best-fit values 

Am^ f M = 2.6 x 1(T 3 eV 2 , sin 2 2^ TM = 1 . (4.29) 

Hence, the best-fit effective atmospheric mixing is maximal. The right panel in Fig. El 
shows the region in the sin 2 2?9ATM~A m ATM plane for — * v T oscillations allowed by 
Super-Kamiokande data |15()j . The left panel in Fig.Qshows the 90% C.L. allowed regions 
for — f v T oscillations obtained in the MACRO and Soudan-2 experiments confronted 
with the corresponding region obtained in the Super-Kamiokande experiment |2S]. The 
right panel in Fig. [7| shows the allowed regions for disappearance obtained in the K2K 
experiment confronted with the allowed regions for — > v T oscillations obtained in the 
Super-Kamiokande experiment |151j . The left panel in Fig. |H1 shows the allowed region 
obtained in Ref. |153j from the combined analysis of Super-Kamiokande atmospheric and 
K2K data. 

Transitions of atmospheric u^s into z/ e 's or sterile states are disfavored. The fraction 
sin 2 £ of atmospheric v^s that transform into sterile v s [y^ — > cos£ v T + sin£ z/ s ) is limited 
at 90% C.L. by [H2] 

sin 2 i< 0.19. (4.30) 

In the next years the MINOS [ZO] experiment will measure with improved precision 
the disappearance of muon neutrinos over a long-baseline of about 730 km. The OPERA 
|155j and ICARUS |156j experiments belonging to the CERN to Gran Sasso program 
(CNGS) [Ell] are aimed at a direct measurement of — > v T oscillation over a similar 
long-baseline of about 730 km. 

4.3 The reactor experiment CHOOZ 

CHOOZ was a long-baseline reactor v e disappearance experiment |157t 1139} 177] which 
did not observe any disappearance of electron neutrinos at a distance of about 1 km from 
the source. In spite of such negative result, the CHOOZ experiment is very important, 
because it shows that the oscillations of electron neutrinos at the atmospheric scale of 
Am 2 are small or zero. This constraint is particularly important in the framework of 
three-neutrino mixing, as will be discussed in Section Therefore, we briefly review the 
results of the CHOOZ experiment. 

The CHOOZ detector consisted in 5 tons of liquid scintillator in which neutrinos were 
revealed through the inverse /3-decay reaction 17 

u e +p-^n + e + , (4.31) 

17 The inverse /3-decay reaction (14.31(1 has been used by all experiments aimed at the detection of 
reactor electron antincutrinos, starting from the Cowan and Reines experiment in 1953 158 , in which 
neutrinos were detected for the first time. The same reaction is used in the KamLAND experiment 
discussed in Section |4~T1 
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Figure 8: Left: Allowed region obtained from the analysis of Super-Kamiokande atmo- 
spheric and K2K data in terms of — > v T oscillations [153J. Right: CHOOZ exclusion 
curves |139j confronted with the Kamiokande allowed regions [81J. 



with a threshold E t ^ = 1.8 MeV. The neutrino energy is measured through the positron 
energy: E = E e + — 1.8 MeV. The detector was located at a distance of about 1 km from 
the Chooz power station, which has two pressurized-water reactors. 

The ratio of observed and expected number of events in the CHOOZ experiment is 

atCHOOZ 

-^ = 1.01 ±0.04, (4.32) 

expected 

showing no indication of any electron antineutrino disappearance. The right panel in 
Fig. El [139 shows the CHOOZ exclusion curves confronted with the Kamiokande allowed 
regions for — > v e transitions [HI]. The area on the right of the exclusion curves is ex- 
cluded. Since the Kamiokande allowed region lies in the excluded area, the disappearance 
of muon neutrinos observed in Kamiokande (and 1MB, Super-Kamiokande, Soudan-2 and 
MACRO) cannot be due to z/ M — > v e transitions. Indeed, — > v e transitions are also 
disfavored by Super-Kamiokande data, which prefer the — > v T channel [154] (there- 
fore, the Super-Kamiokande collaboration did not calculate an allowed region for — > v e 
transitions and the CHOOZ collaboration correctly compared their exclusion curve with 
the regions allowed by the results of the Kamiokande experiment). 

The results of the CHOOZ experiment have been confirmed, albeit with lower accu- 
racy, by the Palo Verde experiment [78] . 



5 Phenomenology of three-neutrino mixing 

The solar and atmospheric evidences of neutrino oscillations are nicely accommodated 
in the minimal framework of three-neutrino mixing, in which the three flavor neutrinos 
i/ e , z/ M , v T are unitary linear combinations of three neutrinos vi y with masses mi, 
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Figure 9: The two three- neutrino schemes allowed by the hierarchy Am§ UN <C Am\ TM . 

7712, m 3, according to Eq. ([2.46)1 . As explained in Section El this scenario is theoretically 
motivated by the see-saw mechanism, which also predicts that massive neutrinos are 
Majorana particles. 

5.1 Three-neutrino mixing schemes 

Figure 01 shows the two three-neutrino schemes allowed by the observed hierarchy of 
squared- mass differences, Awig UN <C Am\ TM , with the massive neutrinos labeled in order 
to have 

AmluN = AT712! , Am^TM ~ lAm^l ~ |Am|j| . (5.1) 

The two schemes in Fig. El are usually called "normal" and "inverted" , because in the 
normal scheme the smallest squared-mass difference is generated by the two lightest 
neutrinos and a natural neutrino mass hierarchy can be realized 18 if mi <C 777.2, whereas in 
the inverted scheme the smallest squared-mass difference is generated by the two heaviest 
neutrinos, which are almost degenerate for any value of the lightest neutrino mass 777-3. 
This is shown in Fig. where we have depicted the allowed ranges (between the dashed 
and dotted lines) for the neutrino masses obtained from the allowed values of Am| UN 
in Eq. ([4.16J) and Am^ TM in Eq. (|4.27j) . as functions of the lightest mass in the normal 
and inverted schemes. The solid lines correspond to the best fit values of Amg UN and 
Am\ TM in Eqs. ()4.18|) and (|4.29j) . respectively. One can see that at least two neutrinos 
have masses larger than about 7 x 10 _3 eV. 

In the case of three-neutrino mixing there are no light sterile neutrinos, in agreement 
with the absence of any indication in favor of active-sterile transitions in both solar and 
atmospheric neutrino experiments. Let us however emphasize that three-neutrino mixing 
cannot explain the indications in favor of short-baseline — > v e transitions observed in 
the LSND experiment [75] . which are presently under investigation in the MiniBooNE 
experiment [119j . 

Let us now discuss the current information on the three-neutrino mixing matrix U. 
In solar neutrino experiments and v T are indistinguishable, because the energy is well 

18 The absolute scale of neutrino masses is not determined by the observation of neutrino oscillations, 
which depend only on the differences of the squares of neutrino masses. 
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Figure 10: Allowed ranges for the neutrino masses as functions of the lightest mass m\ 
and 777,3 i n the normal and inverted three-neutrino scheme, respectively. 



below fj, and r production and v T can be detected only through flavor-blind neutral- 
current interactions. Hence, solar neutrino oscillations, as well as the oscillations in the 
KamLAND experiment, depend only on the absolute value of the elements in the first row 
of the mixing matrix, \U e i\, \U e ?>\ which regulates v e and v e disappearance. Indeed, 

the survival probability of solar electron neutrinos can be written as [159] 

^^(l-l^l 2 ) 2 ^' 2 ^!^! 4 , (5.2) 

where P^^^f^ is the two-neutrino survival probability in matter (J3.50)) calculated with 
the charged- current matter potential Vcc multiplied by (1 — |?7 e 3| 2 ) and d = $i 2 in the 
parameterization ()2.48|) of the mixing matrix. 

The hierarchy Amg UN <C Am\ TM implies that neutrino oscillations generated by 
Am\ TU depend only on the absolute value of the elements in the last column of the 
mixing matrix, |£7 e 3|, I |, \U T s\ because mi and m 2 are indistinguishable. Indeed, 
taking into account also the matter effects in the earth, the evolution equation of the 
neutrino amplitudes is given by Eq. (|3.27j) with Am|j ~ Am^fe, leading to 

l ^ MX) = J2 ( U (* U * P 3 + he he ^Cc] ^a P (x) , (5.3) 

which clearly depends only on the elements U e s, and U t3 of the mixing matrix. 
In order to further demonstrate that only the absolute values |t/ e3 |, |f^3|, \U T $\ are 
relevant, we notice that in the parameterization ()2.48|) of the mixing matrix we have 
Up 3 = \Up 3 \e- iipi35 ^, and Eq. (JE3J) can be written as 

i ^ ^(x) = e-^- (Wl \U P s\ + he he Vcc) e^ 5 <>° ^ ap (x) . (5.4) 

Since the flavor transition probabilities depend on the squared absolute value of the flavor 
amplitudes (see Eq. (|3.26| ). we can change arbitrarily the phases of the flavor amplitudes. 
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Making the change of phase 

Mx) ^ e-^^^ix) , (5.5) 
we obtain the evolution equation 

^ / Awi 2 \ 

i-^ipapix) = z2[ \ U fr\ ' f/p3 ' + S P eS P e Vcc ) ^p(x) , (5.6) 

p ' 

which depends 19 only on \U e3 \, and \U t3 \. 

The only connection between solar and atmospheric oscillations is due to \U^\. There- 
fore, any information on the value of \U^\ is of crucial importance. 

The key experiment for the determination of \U^\ has been the CHOOZ long-baseline 
reactor v e disappearance experiment |157| 1139) 177] , which did not observe any disappear- 
ance at a distance of about 1 km from the reactor source (see Section [Q|) . The negative 
result of the CHOOZ experiment, confirmed by the Palo Verde experiment [7S], implies 
that the oscillations of electron neutrinos at the atmospheric scale are very small or even 
zero. The CHOOZ bound on the effective mixing angle (see Refs. |160) 121"]) 

sin 2 2tf CHO oz = 4 \V eZ \ 2 (l - |f/ e3 | 2 ) (5.7) 

implies that |C/ e 3| is small: 

|f/ e3 | 2 < 5 x HT 2 , (5.8) 

at 99.73% C.L. [161J. Therefore, solar and atmospheric neutrino oscillations are practi- 
cally decoupled and the effective mixing angles in solar, atmospheric and CHOOZ 
experiments can be related to the elements of the three-neutrino mixing matrix by (see 
also Ref. |IE2| ) 



12 

1 - Itfeai 



sin 2 $sun = 1 U f T \ |2 , sin 2 $ ATM = \U^\ 2 , sin 2 $ C hooz = l^l 2 • (5.9) 



Taking into account the best-fit values of "#sun and ^atm in Eqs. (|4.18jl and (|4.29jl . 
respectively, and 

sin 2 4 f HOOZ = , (5.10) 
the best-fit value for the mixing matrix U is 

TT ( 0.84 0.55 0.00 \ . . 

Uu— -°- 39 °- 59 °- 71 • (5.11) 

V 0.39 -0.59 0.71 / v ' 

Using the 99.73% C.L. allowed ranges for ^sun, ^atm and ^chooz given by Eqs. (J4.17)) . 
1)4.28)1 and (j5.8j) . respectively, we have reconstructed the allowed ranges for the elements 
of the mixing matrix: 

. . / 0.76-0.88 0.47-0.62 0.00-0.22 \ 

\U\ ~ 0.09-0.62 0.29-0.79 0.55-0.85 . (5.12) 
1 1 V 0.11-0.62 0.32-0.80 0.51-0.83/ v ' 

Such mixing matrix, with all elements large except £7 e3 , is called "bilarge". It is very 
different from the quark mixing matrix, in which mixing is very small. This difference 

19 A simpler way to obtain the same result is to adopt a parameterization of the mixing matrix in 
which the Dirac phase is associated with the mixing angle $12, which does not contribute to the evolution 
equation l|5.3|l (see Ref. [Ml). 
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Figure 11: Effective neutrino mass nip in Tritium /3-decay experiments as a function of the 
lightest mass ni\ and in the normal and inverted three-neutrino scheme, respectively. 



is an important piece of information for our understanding of the physics beyond the 
Standard Model, which presumably involves some sort of quark-lepton unification. 

An important open problem is the determination of the absolute values of neutrino 
masses. The most sensitive known ways to probe the absolute values of neutrino masses 
are the observation of the end-point part of the electron spectrum in Tritium /?-decay, 
the observation of large-scale structures in the early universe and the search for neutri- 
noless double-/? decay, if neutrinos are Majorana particles (we do not consider here the 
interesting possibility to determine neutrino masses through the observation of supernova 
neutrinos; see Ref. [35] and references therein). 



5.2 Tritium f3- decay 

Up to now, no indication of a neutrino mass has been found in Tritium /9-decay experi- 
ments, leading to the 95% C.L. upper limit |163j 

m /3 <2.2eV (5.13) 

on the effective mass 

obtained in the Mainz |164j and Troitsk [165 experiments. After 2007, the KATRIN 
experiment [166] will explore nip down to about 0.2 — 0.3 eV. Figure ^2 shows the 
allowed range (between the dashed lines) for nip obtained from the 99.73% C.L. allowed 
values of the oscillation parameters in Eqs. ()4.16|) . (|4.17j) . ()4.27j) . ()4.28|) . as a function 
of the lightest mass in the normal and inverted three-neutrino schemes. The solid line 
corresponds to the best fit values of the oscillation parameters in Eqs. f|4.18|) and f)4.29j) . 
One can see that in the normal scheme with a mass hierarchy nip has a value between 
about 3 x 10~ 3 eV and 2 x 10 _2 eV, whereas in the inverted scheme mp is larger than 
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about 3 x 1CT 2 eV. Therefore, if in the future it will be possible to constraint nip to be 
smaller than about 3 x 10~ 2 eV, a normal hierarchy of neutrino masses will be established. 

From Figs. El arid ^2 it is clear that the bound ()5.13j) can be saturated only if the 
three neutrino masses are almost degenerate. In this case the dependence of m k on the 
index k can be neglected in Eq. (j5.14j) . leading to nip ~ m k . Therefore, the upper limit 
for each mass is the same as the one on nip in Eq. (j5.13|) : at 95% C.L. 

m k <2.2eV. (5.15) 



5.3 Cosmological bounds on neutrino masses 

In the early hot universe neutrinos are in equilibrium in the primeval plasma through the 

weak interaction reactions z/z/ ^ e + e~, Ve ^ Ve, v N ^ u N, v e n ^ pe~ , v e p ^ ne + , 
n ^ pe~u e . As the universe expands and cools, the rate of weak interactions decreases. 
When the temperature of the Universe goes below T ~ 1 MeV, the mean neutrino free 
path becomes larger than the horizon 20 and neutrinos practically cease to interact with 
the plasma. At T < 0.5 MeV electron and positron in the plasma annihilate into photons, 
increasing the photon temperature T 7 with respect to the neutrino temperature T v by a 
factor (11/4) 1 / 3 , easily calculated from entropy conservation (see Refs. ^3J|H1|32|)- From 
the well measured temperature of the Cosmic Microwave Background Radiation (CMBR) 
T 7 = 2.725 ± 0.001 K, we infer the neutrino temperature T u = 1.945 ± 0.001 K, and 
kT v = 1.676±0.001 x 10~ 4 eV. As we have seen in Section l5~T| at least two neutrinos have 
masses larger than about 7 x 10~ 3 eV (see Fig.llOjl. Hence, at least two massive neutrinos 
in the present relic neutrino background are non-relativistic. The number density of 
relic non-relativistic neutrinos can be calculated from the Fermi-Dirac distribution (see 

Refs. panning): 

3 C(3) 

nu k ,u k = ^V^, 3 - 0.1827T 3 ~ 112 cm" 3 , (5.16) 

with ((3) ~ 1.202. Their contribution to the present density of the universe (normalized 
to the critical density p c = 3H 2 /87tGn, where H is the Hubble constant and Gn is the 
Newton constant) is given by 

Uh ~ Pc -^94.14eV [bA ° 

where h is the value of the Hubble constant in units of 100 km s -1 Mpc" 1 (the current 
determination of h from a global fit of cosmological data is h = 0.71!°,' 03 t 167j) . The total 
contribution of relic neutrinos to the present density of the universe is given by 168, 1( 



Uvh "94.14eV [bAH) 

It is clear that, just from the need to avoid overdosing the Universe, the sum of neutrino 
masses has to be lighter than about 100 eV. If one further takes h < 0.8 and Q u < 0.1, 



^The horizon is the distance traveled by light from the beginning of the universe. 
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as indicated by astronomical data, one gets a quite restrictive upper bound of about 6 
eV for the sum of neutrino masses. 

An even stronger bound on the sum of neutrino masses follows from more sophisticated 
calculations of structure formation in the early universe. Neutrinos with masses of the 
order of 1 eV or lighter constitute what is called "hot dark matter" , which is dark matter 
that is now non-relativistic, but was relativistic at the time of matter-radiation equality, 
when the contribution of matter and radiation to the density of the universe was equal. 
Since the radiation energy density scales as T 4 and matter energy density scale as T 3 , 
matter started to dominate the density of the universe and structures started to form 
after matter-radiation equality. However, hot dark matter particles did not participate 
to the beginning of structure formation at matter-radiation equality, but streamed freely 
until they become non relativistic. Hence, neutrinos contribute only to the formation of 
structures with size given by the free-streaming distance traveled by neutrinos until they 
become non relativistic. The formation of structures on smaller scales is suppressed with 
respect to a universe without hot dark matter. The absence of such suppression in the 
present astronomical observations of large scale structures (LSS) in the universe allow to 
put a strong upper bound on the sum of neutrino masses (see Refs. |17()| I17l| 1172] ). 

The recent high-precision CMBR data of the WMAP satellite |173| combined with 
the LSS data of the 2dF Galaxy Redshift Survey (2dFGRS) |174j and other astronomical 
data (see Ref. |lb'7| ) allowed the WMAP collaboration to derive the impressive bound 

tt u h 2 < 0.0076, (5.19) 

with 95% confidence, which, using Eq. (|5.18j) . yields 

^m fc <0.71eV. (5.20) 

k 

From the smallness of the squared-mass differences implied by solar and atmospheric 
neutrino data (see Eqs. (|4.1fij) . (|4.27|) and (j5.1jl ). it is clear that the bound (|5.20|) can be 
saturated only if the three neutrino masses are almost degenerate. Therefore, the upper 
limit on each neutrino mass is one third of the bound in Eq. (j5.20j) : 

m fc <0.23eV. (5.21) 

This impressive limit is one order of magnitude more restrictive than the limit ()5.15|) 
obtained in Tritium experiments, reaching already the level of sensitivity of the future 
KATRIN experiment. Let us emphasize, however, that the KATRIN experiment is im- 
portant in order to probe the neutrino masses in a model-independent way. Indeed, 
the cosmological bound relies on several assumptions on the cosmological model and 
some controversial data (see the discussion in Ref. |175j and Ref. |176j for an alternative 
model). Using only the WMAP and 2dFGRS data, the author of Ref. [TTTj found the 
95% confidence limit 

^m fc <2.12eV. (5.22) 

k 

Adding also the Hubble Space Telescope determination of h, the authors of Ref. |178j 
obtained the 95% confidence limit 

^mjfc<l.leV. (5.23) 

k 
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Figure 12: Effective Majorana mass in neutrinoless double-/? decay experiments 

as a function of the lightest mass mi and in the normal and inverted three-neutrino 
scheme, respectively. 



5.4 Neutrinoless double-/? decay 

A very important open problem in neutrino physics is the Dirac or Majorana nature of 
neutrinos. From the theoretical point of view it is expected that neutrinos are Majorana 
particles, with masses generated by the see-saw mechanism (see Section 12. 4j) or by effec- 
tive Lagrangian terms in which heavy degrees of freedom have been integrated out (see 
Section 1231 and Ref. |3E]). 

The best known way to search for Majorana neutrino masses is neutrinoless double-/? 
decay, whose amplitude is proportional to the effective Majorana mass (see Refs. fSl El 

HUSHED]) 



Km) I 



Yl U zk m k 



(5.24) 



The present experimental upper limit on between about 0.3 eV and 1.3 eV, has 

been obtained in the Heidelberg-Moscow |179j and IGEX experiments [180J. The large 
uncertainty is due to the difficulty of calculating the nuclear matrix element in the decay. 
Figure IT2l shows the allowed range for \{m)\ obtained from the 99.73% CL. allowed values 
of the oscillation parameters in Eqs. (|4.1fij) . (j4.17|) . ()4.27j) . ()4.28j) . as a function of the 
lightest mass in the normal and inverted three- neutrino schemes (see also Refs. |18H 1182^ 
183J). If CP is conserved, |(m)| is constrained to lie in the shadowed region. Finding 
\(m)\ in an unshaded strip would signal CP violation. One can see that in the normal 
scheme large cancellations between the three mass contributions are possible and |(m)| 
can be arbitrarily small. On the other hand, the cancellations in the inverted scheme are 
limited, because V\ and 1/2, with which the electron neutrino has large mixing, are almost 
degenerate and much heavier than z/ 3 . Since the solar mixing angle is less than maximal, 
a complete cancellation between the contributions of u% and v 2 is excluded, leading to a 
lower bound of about 7 x IO -3 eV for | (m) | in the inverted scheme (see also Ref. [184j ). If 
in the future I (m) I will be found to be smaller than about 7x 10 -3 eV, it will be established 
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that either neutrinos have a mass hierarchy or they are Dirac particles. Many neutrinoless 
double-/? decay experiments are planned for the future (see Refs. |185| 135 } I187j ) . The 
most sensitive may be able to probe such small values of | (m) \ . 

6 Future prospects 

As reviewed in Section in recent years neutrino oscillations have been established in a 
model-independent way in solar (~ 19a) and atmospheric (~ 7a) neutrino experiments. 
An impressive proof of neutrino oscillations has also been obtained in the KamLAND re- 
actor experiment (~ 4a). As we have discussed in Section^ all these experimental results 
are nicely accommodated in the framework of three-neutrino mixing. Only the LSND 
anomaly [IS] does not fit in this scheme, as already noticed. Its interpretation in terms of 
neutrino oscillation is presently under investigation in the MiniBooNE experiment |119j . 

Putting aside for the moment the controversial LSND anomaly, the immediate fu- 
ture prospects of neutrino oscillation research deal with the precise determination of the 
parameters of the three-neutrino mixing matrix. The actual value of \U e s\, the real un- 
known, will be searched for in the near future with the long baseline accelerator neutrino 
programs MINOS [23 CHE! in the United States, OPERA jM?] and ICARUS dHD! in 
Europe (CERN to Gran Sasso [20]) and JHF |191j in Japan. All these projects use 
conventional neutrino beams. 

However, the evolution of neutrino physics demands new schemes to produce intense, 
collimated and pure neutrino beams. New possibilities have been studied in the last few 
years: neutrino beams from a Neutrino Factory, Beta-Beams and Super-Beams, that here 
we briefly summarize. 

The current Neutrino Factory concept implies the production, collection, and storage 
of muons to produce very intense beams of muon and electron neutrinos with equal 
fluxes through the decays (j4.21j) . Research and development addressing the feasibility 
of a Neutrino Factory are currently in progress. Review studies on the physics reach of 
the Neutrino Factory option are given in Refs. [EUl I193j . Incidentally, we notice that 
anomalous muon decays, due to non standard weak interactions, if responsible of the 
observed LSND effect, could be easily tested at a Neutrino Factory with a short-baseline 
experiment, using a 10 ton detector capable of charge discrimination |125j . 

The Beta-Beam concept, first proposed in Ref. [194j . is based on the acceleration 
and storage of radioactive ions. The /3-decay of these radioactive ions can produce a 
very intense beam of electron neutrinos or antineutrinos with perfectly known energy 
spectrum. The physics reach of a CERN based Beta-Beam and of a Super-Beam+Beta- 
Beam combination is studied in Ref. |195j . In this study a Super-Beam is a very intense 

beam that has a Super Proton Linac (SPL) as injector delivering 10 23 protons on 
target per year with energy of 2.2 GeV (see also Ref. |196j ). The sensitivity to the 
Dirac CP-violating phase (^13 in the neutrino mixing matrix ()2.48|) reachable in a Super- 
Beam+Beta-Beam combination results to be comparable to the sensitivity reachable in 
a Neutrino Factory if sin 2 $13 > 10 -4 . Comparisons between the sensitivities of different 
projects already approved or planned are presented in Ref. |197j . 

A next-generation neutrino oscillation experiment using reactor antineutrinos could 
give important information on the size of the mixing angle $13 |198|l2Ul)j . Reactor exper- 
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iments can give a clean measure of the mixing angle without ambiguities associated with 
the size of the other mixing angles, matter effects, and effects due to CP violation. The 
key question is whether a next-generation experiment can reach the needed sensitivity 
goals to make a measurement of sin 2 2 i i9 13 at the 10 -2 level |1991 200 . 

However, the search for \Uiz\ and CP violation in the lepton sector does not cover 
all the items in today neutrino physics. Let us emphasize that still several fundamental 
characteristics of neutrinos are unknown. Among them, the Dirac or Majorana nature of 
neutrinos, the absolute scale of neutrino masses, the distinction between the normal and 
inverted schemes and the electromagnetic properties of neutrinos. 

In our opinion the most important question in today neutrino physics is: are massive 
neutrinos Majorana particles? This question can be resolved with an affirmative answer 
if neutrinoless beta decay is observed |201) [202 . Many experimental proposals exist 
that will increase dramatically the sensitivity of the neutrinoless double-/? decay search, 
reaching |(m}| ~ 10" 2 eV (see Refs. [THK1 IHKl ITHH1 ITHTj l. 

Besides masses, neutrinos can have magnetic moments (see Refs. [T71 ITB"! IT4*1 I36|). In 
the case of a Dirac neutrino, a large enough magnetic moment could lead to spin preces- 
sion in a transverse magnetic field |203[ I2()4| 12*05. 206j, generating transitions between a 
neutrino with negative helicity and one with positive helicity, which is practically sterile. 
A Majorana neutrino cannot have a magnetic moment, but different Majorana neutri- 
nos can have transition magnetic moments, which could lead to Spin-Flavor Precession 
(SFP). However, SFP is suppressed in vacuum by the mass difference of different neu- 
trinos. In 1988 Akhmedov |207| I208j and Lim and Marciano |209j discovered that the 
mass difference can be compensated by the matter potential for neutrinos propagating in 
a medium, leading to the so-called Resonant Spin-Flavor Precession (RSFP) mechanism, 
which is analogous to the MSW effect. 

The RSFP mechanism was proposed as a possible explanation of the solar neutrino 
problem and was considered a viable possibility for many years. Now we know that the 
solar neutrino problem is due to neutrino oscillations (see Section l4~Tj) . but RSFP could 
still be a subdominant mechanism (see Ref. |21()j ). A possible signature of RSFP would be 
a periodic variability of the solar neutrino flux due to temporal variations of the magnetic 
field or to the inclination of the solar equator with respect to the ecliptic. A recent analysis 
of Super-Kamiokande data put severe limits to possible periodic modulations of the 8 B 
solar neutrino flux |211j (see, however, Ref. |212j for a different point of view). 

If massive neutrinos are Majorana particles, a combination of the RSFP mechanism 
and oscillations can induce besides v e — > u^, u T , u^, v T transitions also v e — > v e transitions 21 
of solar neutrinos |213) I214[ I215j . which can be observed in solar neutrino detectors 
through the inverse /9-decay process (14.31)1 . However, the Super-Kamiokande experiment 
did not find any indication of z/ e 's coming from the sun and established an upper limit 
of 8 x 10~ 3 for the averaged probability of v e — > v e conversion in the energy range 8-20 
MeV (assuming the initial Standard Solar Model 8 B v e flux) [216J. 

In any case, even if at the moment there seems to be no indication of an effect of the 
RSFP mechanism in solar neutrino data, it is important to pursue this line of research, 
because magnetic moments are important properties and the existence of large neutrino 
magnetic moments could give crucial indications on the physics beyond the Standard 

21 Here v e is the conventional name for a neutrino state with almost exact positive helicity, which can 
induce the creation of a positron upon scattering with matter (see the remark at the end of Section \2. 
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Model. 

Direct measurements of neutrino magnetic moments are also planned: a proposal on 
the direct detection of antineutrino-electron scattering with an artificial tritium source 
|217j will search for a neutrino magnetic moment about two orders of magnitude smaller 
than the present-day laboratory upper limit, reaching a sensitivity of about 10 _12 /xb, 
where /zb is the Bohr magneton. 

7 Conclusions 

The recent years have been extraordinarily fruitful for neutrino physics, yielding model- 
independent proofs of solar and atmospheric neutrino oscillations, which have been con- 
firmed, respectively, by the reactor experiment KamLAND and the accelerator experi- 
ment K2K. Taking into account the negative result of the CHOOZ long-baseline reactor 
v e disappearance experiment, the global fit of solar, KamLAND, atmospheric and K2K 
data have provided important information on the neutrino mixing parameters in the 
framework of three-neutrino mixing, which is predicted by the natural versions of the 
see-saw mechanism. 

The only experimental result that cannot be explained in the framework of three- 
neutrino mixing is the controversial indication in favor of short-baseline — > v e tran- 
sitions observed in the LSND experiment. It is very important that the MiniBooNE 
experiment running at Fermilab will check the validity of the LSND indication in the 
near future. If MiniBooNE will obtain a positive result, the investigation of new possibil- 
ities (as four-neutrino mixing, CPT violation, etc.) will become imperative. These new 
phenomena would be very interesting for our understanding of the physics beyond the 
Standard Model. 

Even if the values of some parameters of three-neutrino mixing are determined with a 
precision that was unthinkable a few years ago, still several fundamental characteristics of 
neutrinos remain unknown. Among them the most important are: the Dirac or Majorana 
nature of neutrinos, the absolute scale of neutrino masses, the distinction between the 
normal and inverted schemes, the value of {U^l, the existence of CP violation in the lepton 
sector, the number of light neutrinos and the electromagnetic properties of neutrinos. 
Several existing experiments and future projects are aimed at the exploration of these 
characteristics, which are very important for our understanding of neutrino physics. Their 
determination is likely to shed some light on the new physics beyond the Standard Model. 
Therefore, we think that interesting years lie ahead in neutrino physics research. 
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